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ABSTRACT 


A  simple,  theory  based  upon  traveling  wave  concepts  and  the 
Wiener-Hopf  technique,  is  developed  which  describes  the  current 
distributions  on  tubular,  cylindrical  receiving  and  transmitting  antennas. 

A  close  examination  of  the  conditions  necessary  to  obtain  sufficiently 
accurate  asymptotic  solutions  for  reflected  current  distributions  is 
given  along  with  several  numerical  examples  for  cooboration.  This 
along  with  corresponding  modifications  to  other  relevant  terms  in  the 
traveling  wave  solution  for  a  finite  length  cylindrical  antenna  allow 
for  the  consideration  of  a  much  wider  range  of  cylindrical  antennas 
than  normally  possible  under  the  traditional  thin-wire  appro  ximatio  ns  9  — ~y 

ka  «  l  twd  kh  2  i .  Specific  examples  discussed  include  electrically  ' 
short,  (kh  =  0.4  and  fi(h)  =  2  Jln(2h/a)  =  10),  practical  half-wave, 

(kh  *  ir/2  and  fi(h)  =  2  £n(2h/a)  =  10),  and  electrical  thick,  (ka  =  1 
and  kh  =  3ir),  receiving  and  transmitting  antennas.  Comparisons  with 
existing  theories  in  these  cases  and  others  yield  very  acceptable  agree¬ 
ments.  further,  the  receiving  antenna  formulation  allows  for  an  arbitrary 
angle  of  incidence,  0  <_  0^  <_  ir,  of  the  uniform  plane  wave  and  the 
transmitting  antenna  formulation  gives  excellent  input  conductance  data 
over  an  extremely  wide  range  of  antenna  parameters.  Discussions  are  given 
for  such  related  topics  as  the  current  distributions  on  die  internal  wall 
of  a  cylindrical  antenna,  loaded  cylindrical  antennas  and  the  far  field 
radiation  pattern  of  a  cylindrical  transmitting  antenna. 


I 


Table  of  Contents 


Section 


Page 


1.  Introduction 


2.  The  canonical  Integral  in  cylindrical  antenna  problems  . 


3.  Currents  on  cylindrical  antennas 


3.1  Primary  receiving  current  . 

3.2  Primary  transmitting  current  . 

3.3  Secondary  current  on  a  semi-infinite  receiving  antenna. 

3.4  Secondary  current  on  a  semi-infinite  transmitting 

antenna . 

4.  Numerical  comparisons:  infinite  and  semi-infinite  antennas. 


4.1  Primary  transmitting  current 

4.2  Reflection  coefficient,  RCB^) 

4.3  Reflected  current  distributions 


5.  Approximate  expressions  for  the  external  currents  on  finite 
length  cylindrical  antennas  . 


5.1  Finite  receiving  antenna  . 

5.2  Finite  transmitting  antenna 


6.  Approximate  expressions  for  the  internal  currents  on 
cylindrical  antennas . 


6.1  Internal  current  on  a  semi- inf inite  receiving  antenna  . 

6.2  Internal  current  on  a  semi-infinite  transmitting  antenna 

6.3  Internal  current  on  a  finite  length  receiving  antenna  . 

6.4  Internal  current  on  a  finite  length  transmitting  antenna 

6.5  End  conductance  of  a  finite  length  cylindrical  antenna. 

7.  Numerical  results  for  the  finite  length  cylindrical  antenna.  . 


7.1  Current  distribution  on  a  receiving  antenna  .  . 

7.2  Current  distribution  on  a  transmitting  antenna 


8.  Electrically  short  cylindrical  antennas 


8.1  Short  receiving  antenna  .  . 

8.2  Short  transmitting  antenna 


/  tftiW-  v-tf'--  -Vr>~  f 


9.  Concluding  remarks 


77 


Section 


References 
Appendix  A: 
Appendix  B: 
Appendix  C: 

Appendix  D: 

Appendix  E: 
Appendix  F: 
Appendix  G: 

Appendix  H: 


1 


. 81 

Approximate  solutions  for  the  canonical  integral,  iKO^jz).  84 
Input  conductance  of  an  infinitely  long  cylinder .  90 


The  kernel,  K(a),  and  the  factorized  functions,  K  (a) 

and  K_(a) . +  .  .  .  94 

Limiting  forms  of  the  currents  on  cylindrical  antennas 
near  grazing  incidence  and  near  the  ends . 102 

D.l  Semi-infinite  receiving  antenna  .....  .  102 

D.2  Semi- infinite  transmitting  antenna  .  105 

D.3  Finite  length  receiving  antenna . 106 

D.4  Finite  length  transmitting  antenna  .  Ill 

Approximation  for  the  summation  in  the  expression  for 

the  current  at  the  cylinder  end . 113 

The  end  conductance  of  a  finite  length  cylindrical 
antenna  ...... 1 . . 

Electrically  thick  antennas  on  the  order  of  a  few 
wavelengths  in  length  .  122 

G.l  Electrically  thick  receiving  antenna  .  122 

G.2  Electrically  thick  transmitting  antenna . 124 

Far  field  radiation  from  a  cylindrical  transmitting 
antenna . . 


Accession  For 

'NITS  GRAScI 
CTIC  TAB 

Unannounced 

justificaUc 


Distrihut Ion/  _ 

Availability  Codes_ 

!\vail  and/or 
Special 


1.  Introduction 


As  is  well  known,  thin-wire  conductors  are  commonly  used  as  radia¬ 
tors  in  the  design  of  antenna  systems.  The  radius,  a  ,  of  each  wire 
is  typically  much  smaller  than  its  half  length,  h  ,  which  for  most 
applications  is  of  the  order  of  a  wavelength,  X .  Only  in  limited  situ¬ 
ations,  such  as  the  case  of  probing  an  unknown  field,  will  the  length  be 
much  smaller  than  a  free-space  wavelength,  (i.e.,  2h  «  X  )  ,  or  as  in 
the  case  of  a  trailing  antenna  behind  an  aircraft,  will  the  length  be  much 
greater  than  a  free  space  wavelength,  (i.e.,  2h  »  X  )  .  Consequently 
most  linear  antenna  theories,  both  analytical  and  numerical,  are 
developed  with  an  explicit  or  implicit  assumption  that  a  «  X  and 
2h  >  X/2  ,  which  is  commonly  referred  to  as  the  thin  wire  assumption. 

On  the  other  hand,  theories  not  in  this  general  category,  usually  have 
a  much  more  limited  range  of  application,  such  as  for  the  very  short 
antenna  and  the  very  long  antenna. 

More  recently,  the  time-transient  response,  as  well  as  the  broad¬ 
band  frequency  response  of  a  thin-wire  structure  has  become  a  problem  of 
considerable  importance.  For  instance,  in  order  to  access  the  suscep¬ 
tibility  of  a  long  thin  cylindrical  metallic  enclosure,  one  must  obtain 
statistical  information  concerning  the  performance  of  the  cylinder  as  a 
receiving  antenna,  over  an  extremely  wide  frequency  range  as  well  as  an 
arbitrary  angle  of  incidence  (referring  to  illumination  by  an  incident 
plane  weve) .  Computations  not  only  become  excessive  when  conventional 
theories  are  utilized  because  virtually  thousands  of  responses  are 
needed,  but  also  very  awkward  since  different  methods  have  to  be  used 
in  different  frequency  ranges.  A  similar  statement,  of  course,  can  also 
be  made  for  studying  the  impulse  response  of  an  antenna. 


Beginning  with  Hallen's  [1]  integral  equation  formulation  for 
the  current  on  a  cylindrical  antenna  and  including  the  work  of  many 
others  [2]-  [10],  the  thin-wire  approximations  mentioned  above  have 
nearly  always  been  employed.  Weinstein  [2]  did,  however,  observe  that 
his  final  approximate  solutions,  which  were  derived  under  the  thin  wire 
assumptions,  could  be  applied  to  cylinders  having  larger  values  of  ka 
if  the  electrical  length,  kz,  were  very  much  larger.  In  contrast,  the 
theory  of  King  and  Middleton  [11,  Chap.  II],  however,  which  involves  the 
iterative  solution  of  an  integral  equation  for  the  current  on  a  finite 
length  cylinder,  requires  explicitly  that  the  parameter,  ft  =  2  £n(2h/a), 
to  be  large,  h  being  the  half-length  of  the  cylinder.  Although  the 
parameter,  ft,  relates  only  to  the  physical  length  and  radius  of  the 
antenna,  this  approach  still  requires  the  electrical  radius,  ka  ,  to  be 
small  compared  to  unity  and  the  electrical  length,  kh,  cannot  be  very 
small  nor  very  large.  King  also  developed  a  receiving  theory  [11,  Chap.  IV] 
for  antennas  having  a  large  ft.  A  large  ft  was  also  the  basis  for  two 
electrically  short  (kh  S  1)  antenna  theories  developed  by  King  [11, 

Sec. 11.31  and  IV. 8]  and  [12,  Sec.  3.7]  which  were  developed  by  making 
approximations  relevant  to  the  short  antenna  situation  in  the  integral 
equation  formulation  of  the  problem. 

Another  means  of  analysis  for  the  cylindrical  antenna  problem  is 
the  numerical  method  of  Moments  technique  [13],  which  has  the  capability 
of  computing  antenna  characteristics  without  invoking  the  thin  wire 
approximations.  Realistically,  however,  the  computation  time  is  con¬ 
siderable  if  the  antenna  is  not  thin  or  the  length  is  more  than  a  few 
free  space  wavelengths. 
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In  this  paper,  our  aim  is  to  develop  a  simple  unified  theory  for 
computing  the  broadband  characteristics  of  a  transmitting  and/or  receiving 
antenna  when  the  parameter,  ft  ■  2  £n(2h/a)  is  large.  For  a  typical 
thin-wire  antenna  where  ft  ■  2  !n(2h/a)  *  10  ,  our  theory  is  applicable 
for  antenna  lengths  as  short  as  2h  ~  0.12X  and  as  long  as  2h  =  23A  (where  for 
ft  =  10,  ka  is  almost  equal  to  1),  w'niah  in  terms  of  frequency  covers  well 
over  two  orders  of  magnitude  and  is  more  than  adequate  even  for 
transient  computations.  We  also  show  that  our  formulation  may  be  applied 
to  an  electrically  thick  (up  to  ka  ■  1  )  cylindrical  transmitting 
antenna  or  the  electrically  thick  receiving  antenna  (for  the  angularly 
independent  current)  and  obtain  favorable  agreement  with  existing 
theories  even  when  the  parameter  ft  *  2  £n(2h/a)  is  not  large. 

We  begin  with  a  re-examination  of  the  conditions  necessary  to 
obtain  simple  approximate  solutions  to  cylindrical  antenna  problems  via 
the  Wiener-Hopf  technique.  Section  2  discusses  a  pair  of  canonical 
integrals  which  characterize  cylindrical  antenna  problems.  Approximate 
expressions  for  these  canonical  integrals  are  derived  subject  to  the 
condition,  ft(z)  -  2  2n(2z/a)  >>  | 2n[2kz  sin^  (6.^/2)  ]  j  .  The  angle,  , 
refers  to  the  incident  angle  of  the  incoming  wave  and  is  more  fully 
described  later.  In  Section  3  the  various  currents  on  both  infinite 
and  semi- infinite  cylindrical  receiving  and  transmitting  antennas  are 
given  and  their  relationships  to  the  canonical  integrals  established. 

Data  obtained  from  the  approximate  expressions  is  then  compared  with 
numerically  evaluated  "exact"  data  in  Section  4.  A  most  important 
observation  in  this  section  is  that  the  parameter,  ft(z)  ■  2  £n(2z/a)  , 
in  the  basic  condition  of  our  analysis,  need  not  be  very  much  larger 
than  |2n[2kzsin  (6^/z)]|  ,  especially  when  thicker  antennas  (ka  £  0.1) 
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are  involved.  Utilizing  the  process  of  summing  multiple  reflections, 
approximate  expressions  for  the  receiving  and  transmitting  currents  and 
the  input  admittance  for  finite  length  cylindrical  antennas  are  formu¬ 
lated  in  Section  S.  Expressions  for  the  currents  flowing  on  the  internal 
walls  of  receiving  and  transmitting  tubular  antennas  are  given  in 
Section  6.  In  Section  7,  numerical  results  from  our  theory  for  specific 
antennas  are  compared  with  the  results  of  other  authors  using  different 
approaches,  with  acceptable  agreement  in  all  the  cases  considered.  The 
special  case  of  the  electrically  short  antenna  is  discussed  in  Section  8. 
General  conclusions  as  well  as  extensions  of  our  theory  to  loaded 
antennas  and  the  determination  of  the  far  field  radiation  from  a  trans¬ 
mitting  antenna  are  given  in  Section  9. 

The  exact  integral  expressions  appearing  in  Section  4  are  for  the 
most  part,  based  upon  the  Wiener-Hopf  technique  (see  for  example  Nobel 
[14],  Weinstein  [2]  and  Mlttra  and  Lee  [15]).  The  assumed  time  variation 
is  e  and  the  implied  Fourier  transform  pair  is  given  by. 


and 


F(a) 


f  F(z)  ei0tz  dz 
J  —OO 


-±i: 


F(oOe 


-iaz 


(1) 


F(z) 


da 


(2) 
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.  The  canonical  integral  in  cylindrical  antenna  problems 


As  will  be  shown  later,  the  external  current  distributions  on  both 
the  receiving  and  transmitting  cylindrical  antennas  can  be  written  in  terms 
of  the  canonical  integral. 


U(0.;z)  =  -i  —  (1  -  cos  6 
i  n  3 


(k+a) (k  cos  6^+a)K(a) 


;0<z<ro,  0<6^<tt 


The  contour,  F^,  is  shown  in  Figure  1  and. 


where 


K(a)  =  iirJ0(£a)  HQV  1  <Ca> 


£  ■  /k2  -  a2  ■  i/a2  -  k^ 


6^  is  the  incident  angle  of  the  incoming  current  wave  when  (3)  is  used  to 
describe  a  particular  current  distribution  reflected  from  an  end  of  a  cylinder 
and  z  is  a  numerical  distance  along  the  axis  of  the  antenna,  k  =  2ir/X  and 
n  are  the  plane  wave  wavenumber  and  the  intrinsic  impedance,  respectively, 
of  the  medium  surrounding  the  antenna.  The  antenna  to  be  considered  is 
assumed  to  have  an  infinitely-thin,  perfectly-conducting  wall  concentric  about 
the  z-axis  at  a  radius,  a.  The  suppressed  time  factor  is  exp(-iwt),  where  u 
is  the  operating  frequency  in  radians/sec. 

We  shall  also  find  it  useful  to  define  the  auxiliary  canonical  integral, 
Wtfi^z),  which  is  similar  to  U(61;z)  in  (3)  except  for  the  appearance  in  the 
integrand  of  the  additional  function,  K+(a),  defined  as  the  factor  of  K(a)  in 
(4)  which  is  analytic  and  free  of  zeroes  in  the  upper  half  complex  a-plane. 
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w(ei;z) 


-i  —  (1  -  cos  0  ) 

n  i 


K+(-ot)e 


-iaz 


(k+a)(k  cos  0^+a)  K(a) 


r  da 


;  0  i  z  s  »>,  0  s  <  i 


(6) 


Properties  of  K+(a)  are  discussed  in  Appendix  C.  We  note  that  the 
integral  W(0  ;z)  usually  occurs  in  prohlems  concerning  the  currents  re¬ 
flected  from  the  ends  of  cylindrical  antennas.  As  shown  in  (A8)  of 
Appendix  A,  when  z  >>  a,  we  may  approximate  (6)  by 

W(0.;z)  -  K+(k)U(01;z)  (7) 


where  11(0^  z)  is  our  original  canonical  integral  given  in  (3).  Thus 
with  WtO^z)  given  in  terms  of  U(0^;z),  our  particular  use  of  the  form 
of  W(©i;z)  in  (6)  will  be  limited  to  providing  exact  date  (from  the 
numerical  integration  of  (6))  to  compare  with  the  approximate  solutions 
to  follow. 

Subject  to  the  condition, 


n(z) 


2 


>> 


jfcn[2kz  sin2 


(8) 


an  approximate  solution  to  the  canonical  integral,  U(0i;z),  is  obtained 
to  order  [n(z)]“2  in  Appendix  A.  From  (A17)  of  Appendix  A,  this  approxi¬ 
mate  solution  may  be  stated  as 


U(0i;z)  -  ^  elkZUn(f  (0^ ;  z>  -  iir]  -  £n[f(0i;z)  +  in] } 


(9) 


where 


f(ei;z) 


2Cy  +  y  +  ix/2  +  in(2kz) 


Ei(-lvo> 


(10) 
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is  a  slowly  varying  function  of  z  .  is  defined  as 


C  =  -  in(ka)  -  v  ;  Y  =  0.57721 .  1 

w 


which  is  usually  taken  as  a  large  parameter  in  the  typical  thin  wire  appli¬ 
cation  and. 


f6  •' 

v0<ei,z)  =  2  kz  sin2  ~ 


The  function,  E^,  appearing  in  (10)  is  the  exponential  integral  of  the 
first  kind  defined  in  Equation  5.1.1  of  [16].  We  note  that  the  antenna 
parameter,  ft(z) ,  is  defined  in  the  same  way  as  in  [ 1 1 ]  where  it  has  been 
used  as  a  large  parameter  for  the  iterative  solution  of  the  antenna  prob- 


Another  approximate  form  of  the  canonical  integral,  U(0^;z),  which 
stems  from  a  Taylor  series  expansion  of  (9)  subject  to  the  basic  restriction 
stated  in  (8)  is  given  in  (A18)  of  Appendix  A  and  repeated  here, 


“<V‘>  -f  TtfpTT 


Even  though  (9)  and  (13)  are  equivalent  with  respect  to  the  order  of  approxi¬ 
mation  (i.e.,  [0(z)]-2),  we  shall  find  (13)  to  have  a  more  desirable 

behavior  in  the  near-grazing,  0^  -  0  ,  and  near  the  end,  z  ~  0,  situations. 
Otherwise,  (9)  will  appear  to  be  a  more  accurate  result  than  is  (13)  for 
U(01;z). 

It  is  interesting  to  compare  our  approximate  forms  of  U ( 6 ^ ; z )  to 
similar  expressions  derived  by  other  authors.  For  the  current  on  an  infinitely 
long  transmitting  antenna  (0 -  tt),  Shen,  Wu  and  King  [6]  by  a  semi- 


analytical,  semi-curve  fitting  technique  found  a  result  similar  to  our 
U(6;z)  in  (9),  except  that  the  term,  i,n[kz  +  /(kz)2  +  exp(-2y),  replaces 
our  terms,  i.n(2kz)  +  exp(-i2kz)E^ (-i2kz) .  Thus  for  large  kz  ,  our  approxi¬ 
mate  solution  for  U(ir;z)  in  (9)  and  that  of  Shen,  et  al.  [6,  Eq.  6]  are 
quite  similar.  Weinstein  [2]  found  an  approximate  solution  to  an  integral 
similar  to  (3),  (he  called  it  the  "key"  integral),  but  having  a  different 
coefficient  outside  the  integral.  Apart  from  this  coefficient  (our  approach 
introduces  this  term  at  a  later  time),  Weinstein  obtained,  through  a  com¬ 
plicated  variational  approach,  an  approximate  result  equivalent  to  our 
second  approximate  form  of  U(04;z)  in  (13).  Also,  in  a  more  recent  work 
by  Chang,  Lee  and  Rispin  [17],  a  further  approximation  of  (13)  was  obtained 
and  used  in  a  receiving  antenna  analysis.  However,  the  analyses  of  Shen,et  al. 
Weinstein  and  Chang,  et  al.,  mentioned  above,  all  assumed  the  conventional 
thin  wire  restrictions, 

ka  <<  1  (and  kz  >  1)  (14) 

Although  in  [2],  Weinstein  did  observe,  a  posteriori, that  the  approximate 
form  of  his  "key"  integral  (similar  to  (13))  could  be  used  for  larger 
values  of  ka  if  at  the  same  time,  kz  was  very  much  larger.  Thus,  the 
importance  of  our  work  is  not  so  much  contained  in  the  approximate  formulas 
for  U(9^;z)  in  (9)  and  (13),  but  rather  in  the  realization  of  a  less 
restrictive  condition  (given  in  (8))  for  the  validity  of  these  approximate 
formulas.  In  fact,  it  will  be  shown  in  Section  4  that  the  approximate 
formulas  for  U(6^;z)  in  (9)  and  (13)  yield  remarkably  good  agreement  with 
numerically  obtained  "exact"  results  even  when  ft(z)  is  of  the  same  order 
as  |in(vQ)|  .  Hence,  even  the  "much  greater"  restriction  appearing  in 
(8)  can  be  significantly  relaxed. 
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3.  Currents  on  cylindrical  antennas 

In  this  section,  the  currents  on  Infinite  and  semi- Inf inlte  cylindri¬ 
cal  receiving  and  transmitting  antennas,  given  in  terms  of  the  canonical 
integral,  U(8^;z)  in  (3)  of  Section  2,  are  described. 

3.1  Primary  receiving  current 

The  longitudinal  current  averaged  over  the  circumference  on  an 
infinitely  long  cylindrical  antenna  due  to  a  plane  wave  polarized  in  the 
same  plane  as  the  antenna  and  incident  at  an  angle,  6^,  with  respect  to 
the  cylinder  axis  (which  is  also  the  z-axis  as  shown  in  Figure  2a),  may  be 
written  as  [8,  eq.  10] , 

lj(01,z)  -  E*  V(0^;z)  (15) 


where. 


ikzcos0^ 

e  ;  -  *  <  z  <  00  , 

0  <  6^  <  it  (16) 

Jq  is  the  zero  order  Bessel  function  and  K(a)  is  given  in  (4). 
k  -  2  it/ A  and  n  are  the  plane  wave  propagation  constant  and  Intrinsic 
impedance,  respectively,  of  the  surrounding  medium.  Higher  order  varia¬ 
tions  of  the  z-dlrected  current  with  respect  to  the  azimuthal  angle,  4>, 
and  the  ^-directed  currents  on  the  cylinder  are  not  treated  in  this  report. 
Thus,  while  V(0^;z)  represents  the  total  longitudinal  current  on  an 
infinitely  long  electrically  thin  (ka  <<  1)  antenna  very  well,  it  corresponds 
only  to  the  azlmuthally  uniform  longitudinally  directed  current  on  an 
infinite  cylindrical  antenna  in  general. 


.  ,4.  yfc«»toV 

1  i ,z;  kn  Sin  ei  K(kcos0i) 
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Infinitely  long  tubular  cylindrical  antennas 

a.  receiving;  uniform  plane  wave  incident  at  an  angle,  0^, 
with  respect  tp  the  z  (antenna)  axis. 

b.  transmitting;  delta  function  voltage  source  of  strength 
V.  volts,  located  at  z  ■  z.  . 


In  the  special  case  of  a  thin  cylinder,  i.e.,  ka  <<  1,  the  Bessel 
function  and  kernel,  K(a)  in  (4),  may  be  approximated  by  the  leading 


terms  in  their  respective  small  argument  expressions  to  yield, 

ikz  cos  0 . 


4n 


v<V*)  =  g 


- 

sin  0 

sin  0, 

2C  +  Itt  -  2  In 

1 

i 

w 

2 

► 

y 

(17) 


where  C  is  given  in  (11). 
w 


3.2  Primary  transmitting  current 

The  longitudinal  current  on  an  infinitely  long  hollow  cylinder  due  to 
a  uniform  (with  respect  to  the  azimuthal  angle,  $  )  delta  function  voltage 
source  of  strength,  VQ,  at  z  «  zQ  (see  Figure  2b),  may  be  written  as 
t6,  Eq.  1], 


'iaU-z, 


i-  vn 

n  0 


e 


(k2  -a2)K(o) 


da  ; 


»  <  z  < 


(18) 


where  the  contour,  is  shown  in  Figure  1  and  K(a)  is  given  in  (4). 
Comparing  (18)  with  the  canonical  integral  definition  in  (3),  we  may  write 
the  driven  infinite  cylinder  current  as 

I^ZqSz)  -  VqU(it;  jz-zQ|)  ;  -»  <z<  <*>  (19) 


and  use  either  approximate  form  of  U(n,z)  in  (9)  or  (13)  to  determine 
this  current  provided  the  restriction  in  (8)  is  satisfied.  This  procedure 
however,  does  not  yield  a  good  result  at  the  source  since  the  condition  on 
Q(z)  is  violated.  In  particular,  the  real  part  of  the  current  at  the 
source  needs  to  be  evaluated  very  accurately,  since  physically  it 


corresponds  Co  Che  inpuc  conductance  and  hence  Che  power  chac  can  be  radiaced 
from  Che  ancenna.  To  chis  end,  ic  is  shown  in  Appendix  B,  how  an  approximace 
expression  for  primary  cransmiccing  currenC  similar  Co  chac  of  Shen,  ec  al. 


[6]  may  be  conscrucced.  This  currenC  expression  denoced  as  u8(|*-z0l>  re- 
places  Che  cerm  U(tt  ;  |  z-z^  |  )  in  (19)  and  is  given  by 

ik | 2—2  | 

Us(|z-z0|)  “  “  e  0  {in{fg(  i  z-ZqJ  )  -  iir]  -  £n[fg(|z-z0j  )  +  iir]}  (20) 

where 

fg(|z-zQ|)  =  2CW+  V  +  irr/2  +  £n[  (k|z-Zgi)  +  /(k|z-zQ|)2  +  exp(-2y  -  2g)  ] 

(21) 

C  and  y  are  given  in  (11)  and, 
w 

g  -  33.88  (ka)2  exp  |T  (22) 

From  (18)  and  (20),  che  inpuc  conduccance  of  an  inflnicely  long  cylinder 
is  Chen  given  by 

G  (ka)  =  Re  (U  (0)}  =  Re{[-  in(2C  -  g  -  iir/2)  -  in(2C  -  g  +  i3n/2)]} 

®  8  n  w  w 

(23) 

IC  will  be  shown  lacer  in  SecCion  4  chac  (23)  yields  a  very  good  inpuc  con¬ 
duccance  for  an  inflnicely  long  cylinder  as  Chick  as  ka  ■  1.0.  Also,  we 

noce  that  U  (z)  in  (20)  is  asymptotic  to  both  forms  of  U(ir;z)  in  (9) 

8 

and  (13)  for  large  kz  and  differs  only  in  the  vicinity  of  the  source, 
kz  =  0. 

I 

3.3  Secondary  current  on  a  semi-infinite  receiving  antenna 
The  secondary  current  on  the  external  wall  of  a  semi-infinite  receiving 
cylinder  (see  Figure  3a),  arises  from  the  reflection  of  the  current 


15 


I.  (Sj.z)  in  (15)  from  the  end  of  the  cylinder.  From  a  Wiener-Hopf  analysis 
this  reflected  current  may  be  expressed  as  [17,  Eq.  27], 

Irefl<0i*2>  “  Ee  v<eii°>  §  < 1-cos  01)K+(-kcos  0±) 

K+(-o)  e“iaz 

(k+a)  (kcos  e^ajKOx)  da  ;  0  ~  z  <  (2A) 

ro  0  <  0 .  <  it 

i  “ 

The  contour  is  shown  in  Figure  1,  and  as  previously  noted,  K+(a)  comes 

from  the  factorization  of  K(a)  in  (A)  into  functions  analytic  in  the  upper 
and  lower  halves  of  the  complex  a  plane,  i.e.,  K(a)  -  K+(a)K  (a).  This 
factorization  is  more  fully  discussed  in  Appendix  C.  The  superscript  R  in 
(24)  signifies  the  receiving  situation.  We  may  write  (24)  in  terms  of  the 
auxiliary  canonical  integral  W(8^iz)  in  (6),  and  by  virtue  of  (7),  we  have 
the  approximate  expression, 

Irefl(ei,z)  =  “  E0  v(V°)R(ei)U(6i;z);  0  <  z  <  •  , 

0  <  0±  <  IT  (25) 

where  we  have  defined  the  "reflection  coefficient". 


R(V  "  it  K+(k)K+(-k  cos  01)  (26) 

The  approximate  expression  for  the  reflected  current  in  (25)  is  valid  if  the 
basic  condition  in  (8)  is  satisfied.  We  note  that  the  reflected  current 
distribution  considered  here  is,  as  in  the  primary  receiving  current  dis¬ 
tribution  in  Section  3.1,  the  total  z-directed  current  averaged  over  the 
circumference  of  the  cylinder. 

One  of  the  obstacles,  which  in  the  past  has  prevented  the  practical 


f 


application  of  the  Wiener-Hopf  technique  to  thicker  antennas,  has  been  the 
absence  of  a  tractable  expression  for  K+(a)  in  the  range  -k  <  a  <  k 
which  appears  in  (26).  It  is  shown  in  Appendix  C,  that  a  curve  fitting 
procedure  involving  the  factor,  K+(«),  with  compensation  for  its  dominant 
irregularity,  yields  the  approximate  formula, 


K+<«) 


K°(a) 


|A| Cl  +  Br(|>] 


»  0  <  a  <  k 


K°(a) 


|A| [1  +  Br(J)  +  C(|)  ) 


J  {  -k  <  a  <  0 


where  K+(a)  is  the  small  argument  form  of  K+(a)  based  upon  the  assump¬ 
tions  that  ka  «  1  and  oa  <<  1  [l.  Sec.  38]  given  by 


k%>  -  airr-nr  u  - 

w 


]a|  Is  the  magnitude  of  A  given  by. 


A  -  K^OHiTrJ^kajH^tka)] 


is  given  by 


»r  *  U  G~(kM>  -  R«U2C  +  In  +  in(2)J-l) 


and  is  the  real  part  of  a  more  complicated  function,  B,  given  in  (C19)  in 
Appendix  C.  Here  G^(ka)  is  the  input  conductance  of  an  infinite 
cylindrical  antenna  having  an  electrical  radius,  ka,  for  which  we  have  the 
approximate  formula  given  in  (23).  Appendix  B  gives  a  detailed  discussion 


f 
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of  the  exact  and  approximate  forms  of  G^(ka).  And  finally,  the  coeffi¬ 
cient  C  is  given  in  terms  of  | A |  and  by. 


1  -  1 A  | 2  ( 1  -  Br2  ) 
j  A  |  2  (1  +  Br) 


(31) 


Although  (27)  is  basically  a  curve- fit  solution  for  K+(a)  in  the  range, 
-k  <  a  <  k,  the  coefficients  j  A |  and  were  obtained  in  much  the  same 

manner  as  those  in  a  two-term  Taylor  series  expanion  of  K+(ot)/K+(a)  in 
the  upper-half  of  the  complex  ot-plane.  The  coefficient,  C,  was  obtained 
by  requiring  that  the  approximate  constructed  quantity  K(ct)  =  K+(a)K  (a) 
(see  (C5)  and  (C6)  of  Appendix  C)  using  (27)  have  the  same  limiting  form 
as  the  exact  K(a)  in  (4)  as  a  t  k. 

3.4  Secondary  current  on  a  semi- inf inite  transmitting  antenna 
The  secondary  current  eminating  from  the  end  at  z  *  0  of  a  semi¬ 
infinite,  0  <  z  <  “>,  cylinder  having  a  delta  function  voltage  source  of 
strength  at  z  *  z^  is  usually  approximated  by  the  reflection  of  a 

wave  incident  at  0^  =  it  [6],  as  illustrated  in  Figure  3b.  Hence,  from 
(24)  we  may  write, 

Irefl^Z0’Z^  “  -I^(z0;0)R(tt)U(ti;z)  ,  0  <  z  <  ■»  (32) 

where  we  have  replaced  the  Receiving  incident  current,  E*  V(0^;O)  by 

T 

the  transmitting  incident  current,  Iw(z  ;0). 
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4.  Numerical  comparisons:  infinite  and  semi-infinite  antennas 

4.1  Primary  transmitting  current 

As  discussed  in  Section  3.1,  we  shal L  use  the  modified  Shen,  et  al. 

[6]  formula  denoted  as  Us(z)  in  (20)  for  the  primary  current  on  a 
cylindrical  transmitting  antenna.  And  since  Shen  has  already  compared  his 
approximate  expression  with  numerically  "exact"  data  in  [6]  for  values  of 
ka  up  to  0.08  with  good  agreements,  we  shall  only  consider  cases  in  which 
0. 1  <  ka  <  1.0  to  justify  the  extension  of  the  theory  to  this  range.  In 
Figure  4,  we  show  the  real  and  imaginary  components  of  the  current  distri¬ 
bution  on  an  infinitely  long  cylinder  as  predicted  by  the  modified  Shen 
formula  in  (20)  with  V  “  1  volt,  for  the  particular  values  of  the  electri¬ 
cal  radii,  ka  ■  0.1,  0.5,  and  1.0.  "Exact"  data  for  these  cases  obtained 
from  the  numerical  integration  of  (18)  is  also  shown  in  Figure  4  (as 
circles).  And  it  may  be  observed  that  the  real  component  of  the  current 
distribution  predicted  by  (20)  compares  very  favorably  with  the  exact 
numerical  data  over  the  entire  range  of  kz  shown  especially  for  the 
smaller  values  of  ka.  The  imaginary  component  of  the  current  distribution 
predicted  by  (20)  compares  favorably  with  the  exact  data  only  when  the 
ratio,  2z/a,  somewhat  exceeds  unity. 

As  mentioned  earlier,  the  purpose  of  Shen's  and  our  curve-fitting  pro¬ 
cedures  leading  to  (20)  for  the  primary  transmitting  current  was  to  obtain  a 
good  value  for  the  real  part  of  the  current  at  the  source,  i.e.,  the  input 
conductance  of  an  infinitely  long  cylindrical  antenna.  To  demonstrate  the 
level  of  success  attained  in  this  respect  we  offer  Figure  5,  which  shows 
the  input  conductance  of  an  infinitely  long  cylinder  as  obtained  from  the 
real  part  of  Ug(0)  in  (23)  and  the  "exact"  numerically  evaluated  input 
conductance  from  the  exact  integral  expression  stated  in  (Bl)  of  Appendix 


PRIMARY  TRANSMITTING  ANTENNA  CURRENT 
IM AG (ma/V )  REAL (mo/V) 
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Figure  4.  Current  distribution  on  an  infinitely  long  cylindrical  trans¬ 
mitting  antenna  with  a  delta  function  voltage  source  at  z*Zq“0. 

_ _  Approximate  distribution  from  the  modified  Shen  [6] 

formula  denoted  by  U  (|z-Zq|)  in  eq.  (20). 

O  "Exact"  numerically  evaluated  data  from  -*q.(18). 

Note:  The  traveling  wave  factor,  e1^2,  has  been  su,  ..res  .id. 
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B  as  functions  of  ka  over  the  range  10-4  <  ka  <  1.  Obviously,  excellent 
agreement  is  obtained.  In  fact,  the  error,  which  is  also  shown  in  Figure 

5,  never  exceeds  2%  over  the  entire  range. 

4.2  Reflection  coefficient  R(0.) 

i 

The  behavior  of  the  "exact1'  numerically  evaluated  (using  the  formula 
of  Mittra  and  Lee  (15,  Sec.  5-2. (3)])  K+(a)  is  shown  in  Figure  6  as  a 

function  of  a  in  the  range  -k  <  a  <  k  for  the  specific  cases  ka  = 

0.01,  0.05,  0.1,  0.5,  and  1.0.  This  variation  in  a  when  a  =  -k  cos  6 
corresponds  to  the  range  0<  8  <  u.  The  behavior  of  our  approximate  form 
of  K+(a)  in  (27)  is  so  close  to  the  exact  we  have  not  included  this  data 
in  Figure  6  but  have  elected  to  show,  in  Figure  7,  the  error  between  the 
approximate  and  "exact"  values  of  K+(ot)  for  the  same  range  and  set  of 
parameters  as  those  in  Figure  6.  The  magnitude  and  phase  error  illustrated 
in  Figure  7  is  seen  to  be  quite  small,  typically  below  1%  and  ±5°,  respec¬ 
tively.  And  it  should  be  noted,  that  this  magnitude  error  is  many  times 
smaller  than  the  magnitude  error  of  the  normally  accepted  small  argument 
approximation,  K^(a)  in  (28).  For  example,  at  ka  =  0.01  K^(k)  differs 

from  the  exact  value  of  K+(k)  by  about  1.5%,  while  our  approximate  form 
of  K+(k)  from  (27)  possesses  an  error  of  less  than  0.1%.  And  as  the 
value  of  ka  increases,  the  error  in  K+(k)  increases  quite  rapidly, 
reaching  over  200%  at  ka  ■  1. 

Obviously,  the  quantity  of  more  crucial  importance  is  the  so-called 
"reflection  coefficient",  R(8^)  in  (26).  Figure  8  shows  the  magnitude 
and  phase  of  R(0^)  as  calculated  using  the  "exact"  numerically  determined 
values  of  K+(k)  and  K+(-k  cos  0^)  (again  from  the  formula  of  Mittra 
and  Lee  [15,  Sec.  5-2. (3)])  as  a  function  of  ka  over  the  range 
10-4  <  ka  <  1  for  the  incident  angles  8^  *  ti/36,  u/4,  ti / 2 ,  3ir/4  and  it. 


/ 


PHASE  (deg)  MAGNITUDE 


Figure  8.  Magnitude  and  phase  of  the  "reflection  coefficient," 

R^)  *  (n/2TC)K+(k)K+(-k  cos  0^) ,  calculated  from  the  exact 
formula  for  K+(a)  in  [15,  Sec  5-2. (3)]  (See  eq.  (Cll)  of 
Appendix  C  of  this  report.) 


11(6^  ”  (n/2ir)K+(k)K+(-k  cos  0^),  calculated  using  the 


approximate  formula  for 


K+(a)  in  (27). 
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Data  for  R(6i>  using  the  approximate  formula  for  K+(a)  in  (27)  is  not 
included  in  Figure  8  because  of  the  very  close  agreement  it  has  with  the 
exact  data.  Instead,  ve  show  the  error  of  this  approximation  with  regard 
to  the  exact  in  Figure  9  for  the  same  range  and  set  of  parameters  as  in 
Figure  8.  The  magnitude  error  is  seen  to  be  at  most  about  3%  and  typically 
much  less  while  the  very  small  phase  error  is  never  more  than  ±5°. 

4.3  Reflected  current  distributions 

Denoting  the  reflected  current  due  to  a  unit  incident  current  of  the 

N 

form  exptikz  cos  9.]  as  I  ,,  (0.;z)  ,  we  have  from  (6)  and  (24)  the  ex- 
l  reti  i 

pression. 


Irefl(8i*z)  " '  it  M"k  cos  V  W(VZ)  ;  0  <  *  £  - 

0  <  0i  <  n  (33) 

in  both  the  receiving  and  transmitting  situations.  We  note  that  (33)  is 
an  exact  expression  for  the  normalized  reflected  current  in  the  receiving 
situation  (0  <  0^  <  it)  and  is  a  very  good  approximation  for  the  normalized 
reflected  current  in  the  transmitting  situation  (O^tt)  when  the  delta 
function  voltage  source  is  located  sufficiently  away  from  the  end.  From 
(7)  and  (26),  the  approximate  form  of  (33)  is  given  by 

I^efl(®i«z)  *  -R(01)U(6i;z)  ;  0  <  z  s  ” 

0  <  0  <  *  (34) 

To  demonstrate  the  accuracy  attained  with  our  approximate  formulas. 

Figures  10-14  show  the  behaviors  of  the  ’'exact”  reflected  currents  in  (33) 
(with  K+(-k  cos  6^)  numerically  determined  using  the  formula  of  Mittra 
and  Lee  [15,  Sec.  5-2.13]  and  W(0^;z)  in  (6)  numerically  integrated)  and 
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Legend  for  Figures  1 0— 1 4 


O  6^" 

O  -n/2 
O  -»/4 

£.  =ti/ 36 


N 

"Exact"  numerically  determined  I^^^e^jz)  from  (33) 
found  by  using  the  formula  of  Mittra  and  Lee  [15,  Sec. 
5-2. (3)]  for  K+(a)  and  by  the  numerical  integration  of 
W(0^;z)  given  in  (6). 


N 

Approximate  form  of  I  f  (h.;z)  given  in  (35),  with 
R(61)  from  (26)  determined  by  the  approximate  K+(a) 
formula  in  (27)  and  the  approximate  formula  (13)  used  for 
U(6i;z). 


N 

Approximate  form  of  given  in  (35),  with 

R(Q^)  from  (26)  determined  by  the  approximate  K+(a) 
formula  in  (27)  and  the  approximate  formula  in  (9)  used 
for  U(01;z). 


l.k.z 

Note:  In  Figures  10-14,  the  traveling  wave  phase  factor,  e  , 
has  been  suppressed  to  aid  in  improving  the  clarity  of 
the  Information  presented. 
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Figure  10.  Magnitude  and  phase  of  the  current  reflected  from  the  end 
of  a  semi-infinite  tubular  cylinder  where  ka  *  0.01. 

(See  accompanying  Legend  for  further  details.) 
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Figure  14.  Magnitude  and  phase  of  the  current  reflected  from  the  end 
of  a  semi-infinite  tubular  cylinder  where  ka  ■  1.0.  (See 
accompanying  Legend  for  further  details.) 
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the  approximate  reflected  currents  in  (3'.)  using  (27)  for  K+(a)  in  (26) 
for  R(0i)  and  using  both  (9)  and  (13)  for  iKe^jz)  as  function  of  kz 
in  which  0^  =  n/36,  n/4,  tt/2,  and  n  for  the  values  of  ka  equal  to 
0.01,  0.05,  0.1,  0.5,  and  1.0,  respectively.  We  note  that  in  every  case, 
the  data  obtained  from  the  use  of  (9)  for  U(6^;z)  appears  to  be  closer 
to  the  numerical  data  than  does  results  using  (13)  for  U(0  ;z).  This  is 
somewhat  misleading,  since  in  the  finite  length  cylinder  situation  where 
multiple  reflections  of  currents  from  the  end  are  characterized  by  waves 
incident  at  an  angle,  0^  = n  ,  and  subsequently  summed  (see  Section  5) 
slightly  better  results  are  obtained  witu  (13)  used  for  U(0.;z).  This 
apparent  incongruity  must  be  a  result  of  the  summation  procedure  producing 
an  error  which  is  more  compensitive  for  the  error  in  (13)  for  U(@i;z) 
than  it  is  for  the  error  in  (9).  A  more  detailed  clarification  of  this 


point  will  be  forthcoming. 
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Section  5.  Approximate  expressions  for  me  external  currents  on  finite- 
length  cylindrical  antennas 

Expressions  for  the  currents  on  finite  length  cylindrical  antennas 
are  constructed  by  summing  the  primary  and  subsequent  secondary  currents 
reflected  from  the  ends  of  the  antenna. 

5.1  Finite  receiving  antenna 

Our  theory  can  now  be  applied  to  the  finite  length  receiving  antenna 
with  the  understanding  that  only  the  average  (over  the  circumference)  z- 
directed  current  is  obtained.  As  noted  by  Kao  [ 1 8 ]  specifically  for  normal 
incidence  of  the  plane  wave,  this  zero-order  current  is  not  coupled  to  any 
higher  order  variations  of  the  current  with  respect  to  the  azimuthal  angle, 
<p  ,  and  may  be  considered  independently  from  these  higher  order  currents. 
Rispin  and  Chang  [ 1 9 J  have  also  noted  this  to  be  true  for  arbitrary  polari¬ 
zation  and  arbitrary  incidence  of  the  uniform  plane  wave. 

The  constituative  currents  on  a  finite  length  (~h  <  z  s  h),  cylindri¬ 
cal  receiving  antenna  with  radius,  a,  are  shown  pictorially  in  Figure  15. 
Beginning  with  the  plane  wave  induced  primary  current,  E^V(0.;z),  shown  in 

U  1 

Figure  15a,  the  reflections  of  this  current  from  the  end  at  z  =  -h  and 
the  z  =  +h  end  are  determined  to  be  -V  (n-0  jhlRCO^)!!  (G  ,h+z)  and 
-V(0  ;h)R(n-ei)U(u-ei>h-z),  respectively,  as  illustrated  in  Figure  15b. 
These  reflected  currents  then  propagate  toward  opposite  ends  of  the  cylin¬ 
der  (analogous  to  waves  incident  at  an  angle  it  with  respect  to  a  parti¬ 
cular  end)  at  which  point  they  reflect  again  as  — V (tr— 8  ;h)R(0  )U(0^  ;2h) 
R(n)U(0;h-z)  and  -V(0^  ;h)R(n— 0^)L'(tt  —  0^  ;  2h )  R  (ti  )  U  (tj  ;h+z)  ,  respectively. 
Continuing  this  procedure  leads  to  an  inf inite  number  of  reflected  currents 
eminating  from  each  end  of  the  cylinder,  as  suggested  in  Figure  15c.  The 
infinite  series  expressing  the  current  reflected  from  a  particular  end  of 
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the  cylinder  are  in  the  form  of  simple  geometric  series  which  may  be  readily 
summed.  Hence,  we  arrive  at  the  following  expression  for  the  total  exter¬ 
nal  current  on  the  finite  length  cylindrical  receiving  antenna, 

IR(0.;z)  =  E* { V (0 . ; z)  -  V(n-0  .  ,h)R(  .  )U(6.;h+z)  -  CR(u-e  .  )R(tt)U(m  ,h+z) 

1  D  1  1  L  1  1 

-  V(e.;h)R(7i-ei)U(n-0i,h-z)  -  CR(0,.)R(ti)U(t!  ,h-z) }  (35) 


where 


CR(0.) 


[V(0  ;h)R(*-0 .)U(x-0.  ;2h)R(n)U(ir;2h)-V(ir-0  ;h)R(C  )U (0 ± ; 2h) ] 


1  -  [R(t  )U  (it  ;  2h)  ] : 


(36) 


represents  the  total  incident  current  (with  an  analogous  wave  incidence  of 

0^  =  n)  upon  the  end  z  =  +h  due  to  current  reflections  eminating  from  the 

end  at  z  =  -h.  C  (ir-0  )  has  a  similar  interpretation  with  the  ends 

interchanged.  The  terms  involving  R(6  lU^.h+z)  and  ;h-z) 

represent  the  initial  reflections  of  the  primary  current  wave  incident  at 

the  angles,  0^  and  it  -  0^,,  respectively.  Thus,  except  for  the  primary 

term,  V(0^;z),  all  the  other  terms  in  (35)  represent  reflected  currents 

from  the  ends  of  the  cylinder.  Our  expression  for  the  receiving  antenna 

current  in  (35)  agrees  in  form  with  that  of  Weinstein  [2]  and  can  be  shown 

to  be  consistent  with  our  earlier  result  in  [17]  under  the  conventional 

thin  wire  approximations.  A  complete  formal  agreement  between  our  result 

and  that  of  Shen  [7]  occurs  only  when  the  terms  U(0^;z)  and  U(ir-6^,z) 

in  (35)  are  approximated  by  U  (z)  in  (20)  with  the  constant,  g,  deleted. 

s 

The  approximation  of  these  terms  in  this  manner  is  implicit  in  Shen's  [7] 


analysis. 


The  limiting  form  of  the  current  on  .1  finite  length  receiving  antenna 
as  the  angle,  6^,  approaches  grazing  incidence,  i.e.,  ft  -+■  0  or  »  ,  based 
upon  both  approximate  forms  of  U(0  ;z)  Ln  (9)  and  (13)  is  discussed  in 
Appendix  D.  And  it  is  found  that,  while  our  theory  is  not  expected  to  be 
valid  in  this  range  because  of  the  apparent  violation  of  the  restriction, 
ft(z)  >>  |ln(vQ)|  in  (8),  the  approximate  form  of  U(0^;z)  in  (13)  actually 
produces  the  very  physically  acceptable  result  of  a  vanishing  current  as 
0  -K)  or  it.  Also  a  smaller,  magnitude-wise,  result  for  the  current  near 
the  ends  of  a  cylinder  for  a  fixed  incident  angle,  0  is  obtained  in 
Appendix  D,  when  (13)  is  used  for  U  ( 0 ^ ;  ii  3  rather  than  (9).  These  con¬ 
siderations  are  very  important  in  the  cases  when  the  incident  angle  K 
is  near  grazing,  i.e.,  6^  =  0  or  it,  and  when  the  length  of  the  antenna 

becomes  electrically  short. 

And,  finally,  we  note  the  symmetrical  behavior  of  (35)  with  respect 
to  the  incident  angle  of  the  uniform  plane  wave  and  the  position,  z, 

IR(0.;-z)  =  IR(tt-6.;+z)  (37) 

5.2  Finite  transmitting  antenna 

In  much  the  same  manner,  the  current  on  a  finite  length  (-h  <  z  <  +h) 
cylindrical  transmitting  antBnna  of  radius,  a,  due  to  a  delta  function 
voltage  source  of  strength,  Vq,  at  z  =  z^  (see  Figure  16)  may  be 
expressed  in  terms  of  a  primary  current  and  the  multiply  reflected  currents 
from  the  ends.  Figure  16a  illustrates  the  primary  current,  which  we  shall 
approximate  by  U  (|z-Zq|)  from  (20),  eminating  from  the  delta  function 
voltage  source  at  z  =  z^.  These  waves  are  incident  upon  the  ends  of  the 
cylinder  at  an  angle  of  tt  respective  to  the  particular  end.  Hence,  the 
initial  reflections  of  the  primary  current  from  z  =  -h  and  z  =  +h  are. 
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-V„U  (tt,  h+z„)R(n  )U  (tt  ,h+z)  and  -V.U  (tt  ,  h-z„)R(-n)U(n  ,h-z)  ,  respective]  v  ,  as  is 
0  s  0  0  s  0 

shown  in  Figure  16b.  The  reflections  of  these  currents  from  the  respec¬ 
tive  opposite  ends  and  the  subsequent  reflections  which  follow  (Figure 
16c)  lead  to  a  pair  of  infinite  series,  which  are  again  summable.  The 
final  result  for  the  transmitting  current  distribution  is  given  by. 


1  (Z0;Z)  =  V0{Us(iz_Z0|) 


Ug(h+Zg)R(ir)U(Ti  ;h+z)  -  C  (h+Zg)R(n  ) U (tt  ;h+z) 


Ug  (h-zQ)R (ir)U  (it  ; h-z )  -  C  (h-z  )R(n  )U  (tt  ;h-z)  }  (38) 


where 


T 

C  (z) 


iMz)  [R(TT)U(Tr;2h)]2  -  II  (2h-z)R(Ti)U  (it  ;  2h) 
1  -  [R(tt)1'  (tt  ;  2h)  ] 2 


represents  the  sum  of  the  currents  incident  upon  the  z  =  -h  and  +h  ends 

of  the  antenna  when  z  is  taken  as  h+z  and  h-z  ,  respectively,  due  to 

o  o 

current  reflections  eminating  from  the  opposite  end.  Note  that  the  initial 
reflection  of  the  primary  current  from  each  end  is  explicitly  stated  in 
(38),  the  overall  form  of  the  transmitting  current  expression  being  the 
same  as  that  for  the  receiving  current  in  (35).  Our  transmitting  current 
expression  in  (38)  can  be  shown  to  be  equivalent  in  form  to  those  of  many 
other  authors  [l.  Sec.  35.7],  [2],  [6]  and  others. 

However,  unlike  the  expressions  of  these  authors,  our  expression  is 
more  general  and  flexible,  since  wi  claim  it  may  be  used  for  electrically 
short  as  well  as  electrically  thick  antennas  as  long  as  the  basic  condi¬ 
tion  in  (8),  O(z')  =  2  ln(2z'/a)  <■-'  lln(2kz')|  is  satisfied  (note  here 


Illustration  of  the  multiple  reflection  concept  as  applied  to  the  finite 
length  cylindrical  transmitting  antenna. 

a.  Primary  transmitting  current  distribution  from  (19)  and  (20). 

b.  Initial  reflection  of  the  primary  current  from  each  end  from  eq.  (32) 

c.  Subsequent  multiple  reflections. 


40 


z’  refers  to  the  distance  to  the  source,  z-z^,  and  the  distances  to  the 
cylinder  ends,  h+z  and  h-z)  and  appropriately  accurate  values  of  R(tt) 
are  used.  And  in  this  report,  we  often  Lake  the  nominal  measure,  z'  -  h. 

We  note  the  symmetry  in  (38)  with  respect  to  the  source  and  observa¬ 
tion  points ,  i. e. , 

IT(zq;-z)  =  IT(-zq;z)  (40) 

An  approximate  formula  for  the  input  admittance  of  an  asymetrically 
driven  cylindrical  antenna  of  length,  2h,  obtained  by  setting  z=zQ  and 
Vq  *  1  volt  in  (38)  is  given  by 

Y  -  G  -  iB  *  U  (0) 
in  s 

-  Ug(h+z0)R(Tr)U(ir;h+z0)  -  CT(h+z0)R(ir)U(ir ;h+zQ) 

-  Us(h-z0)R(w)U(*;h-z0)  -  CT(h-z0)R(ir)U(n;h-z0)  (41) 

Here  G  is  the  input  conductance  and  B  is  a  "relative"  input  suscep- 
tance.  The  qualification  to  a  "relative"  input  susceptance  is  necessary, 
due  to  the  fact  we  employ  a  delta  function  voltage  source  for  the  excita¬ 
tion  and  the  mathematically  predicted  behavior  of  the  imaginary  part  of 
the  input  current  for  this  excitation  should  exhibit  a  logarithmic  singu¬ 
larity  [ 20 ]  and  [21],  This  singularity  would  indicate  an  infinite  capaci¬ 
tance,  the  so-called  "knife-edge  capacitance"  [22].  However,  the  parti¬ 
cular  way  in  which  the  primary  current  term,  Ug(z),  was  derived  (discussed 
in  Section  3.2)  doe3  not  allow  the  possibility  of  such  a  singularity  in 
this  current  at  the  source.  In  a  realistic  sense,  though,  this  slice 
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capacitance  is  due  to  an  idealization  in  the  mathematical  model  rather 
than  a  physically  occurring  phenomenon  in  the  practical  situation  and  in 
general  does  not  pose  any  difficulties  in  experimental  studies.  Thus, 
the  absence  of  such  a  singularity  in  our  formulation  is  not  unwelcomed. 
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6 .  Approximate  expressions  for  the  internal  currents  on  cylindrical 
antennas 

Thus  far  our  theory  has  considered  only  the  external  current  distri¬ 
butions  on  cylindrical  antennas,  hence,  it  is  appropriate  at  this  time  to 
include  a  complementary  discussion  of  the  current  distributions  on  the 
internal  walls  of  receiving  and  transmitting  cylindrical  antennas.  By 
combining  the  external  and  internal  current  distributions,  the  total  cur¬ 
rent  on  the  antenna  may  be  found.  But  a  more  important  use  of  a  knowledge 
of  the  internal  current  occurs  in  some  electromagnetic  compatibility 
studies  where  it  is  desirable  to  know  the  amount  of  penetration  into  a 
long  thin  metallic  enclosure.  In  many  cases,  the  penetration  is  into  the 
end  of  a  cylinder  and  one  needs  to  know  the  induced  current  on  the  inter¬ 
nal  wall  of  the  cylinder. 

6.1  Internal  current  on  a  semi-infinite  receiving  antenna 

The  TM_  mode  currents  on  the  internal  wall  of  a  semi-infinite 
un 

(0  <  z  <  “)  cylinder  due  to  a  plane  wave  incident  at  an  angle,  0^,  are 

easily  determined  by  a  Wiener-Hopf  analysis  [17,  Eq.  27]  to  be  given  by 

(I*  (6  ;*)}  =  eJ  V(0  ;0)  (1-cos  . )K  <-k  cos  0  )] 

S°°  I  Int  0  1  JlJ  l  +  i 

K+(-a)e-ioz 

j.  (k+a)(kcoset+u)K(a)  d“  ;  0  "  z  "  °° 

1  0  <  0  <  n  (42) 

where  E*V(0.;O)  is  the  incident  current  at  the  end  and  is  given  in  (15). 

0  I 

The  contour,  Tj,  is  shown  in  Figure  1  and  K+  is  the  "plus"  factor  of 
K(a)  in  (4),  which  is  discussed  in  Appendix  C.  Since  the  contour, 


encloses  only  simple  poles  of  [K(a)]_1,  the  integral  may  be  easily 


evaluated  and  (42)  may  be  written  in  the  form 


us.(V2),i„t-EJv(V0>  \  T0„(V°  Y°“Z 

n=*l 


(43) 


where  T  is  a  transmission  coefficient  given  by. 
On 


u  (Yn n“ik) 

T  (0.)  =  +i  —  ( l-cos0  ,  )K  (-kcos0  . )  - - — - — — r-  K(iy  )  (44) 

On'  i  2  l  +'  V  yn  (v  +ikcos6 . )  +'  Y0n  ' ’ 

On  On  l 


and. 


iy0n  =  i  /(P0n/a)2  "  k2  <«) 

is  the  propagation  constant  of  the  TM^  circular  waveguide  mode.  And 
finally,  pQn  is  the  n— ^  ordered  zero  of  the  Bessel  function,  JQ. 

Several  approximations  are  possible  to  allow  us  to  state  the  internal 
current  in  a  more  convenient  form.  The  first  of  which  is  from  the  approxi 
mate  splitting  of  the  asymptotic  form  of  the  kernel,  K(a),  for  large  aa 
and  is  given  by 

K+(ot)  *  /i/ (k+a)a  ;  for  aa  large  (46) 

Numerical  data  comparing  (46)  with  the  exact  value  of  K+(a)  from  the 

formula  of  Mittra  and  Lee  [ 1 5 ,  Sec.  5-2. (3)]  has  shown  good  agreement 

— y  2 

for  a  >  iy^j  up  to  ka  *  1.0  .  Also  for  e  n  <<  1,  the  infinite  sum 
may  be  truncated  at  n=N  and  the  subsequent  loss  of  information  for  the 
smaller  values  of  z  may  be  somewhat  compensated  for,  by  approximating 
the  summation  in  (43)  at  z-0  using  relevant  Taylor  series  expansions 
for  ka  s  1  in  the  manner  described  in  Appendix  E.  The  summation  in  (43) 
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may  then  be  approximately  written  as 


Ji  T°" ' 


-T0nz 


2  O-cosG  )K  (-kcose  )  £  S  (e  )  (ika)m;  z=0 

ra=  1 


I  T0n(lll)e 


— Y  z 
On 

;  e  <<  1 


where 


S  (6  )  =  0.5831 


S2(01)  =  -0. 1364  +  cos  6.] 
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S3(e.)  =  -0.0498  [-£  -  j  cos  0.  -  cos2e .  ] 


<e i )  =  °*0198  +  ^  cos  0.  -  1  cos20i  -  cos3ei] 


which  is  sufficiently  accurate  for  most  engineering  applications  up  to 
ka  *  1. 


6.2  Internal  currents  on  a  semi-in!  Lnite  transmitting  antenna 

The  current  which  penetrates  into  the  end  of  a  semi-infinite 

(0  <  z  <  00 )  cylindrical  transmitting  antenna  having  a  delta  function 

voltage  source  of  strength,  V  volts,  at  z=z„  is  associated  with  TM 

0  On 

circular  waveguide  modes  and  may  be  written  in  an  analogous  manner  with 
respect  to  the  receiving  case  as, 


{IL(z>}<„,  “  II<zn5°) 

s 00  int  00  0 


~  K  <k) 

7T  + 


v  /  v  -iuz 

K+(-a)e 

(k:’-a2)K(a) 


da  ;  0  <  z  <  »  (52) 
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where  I^ZgjO)  is  the  incident  current  from  (19)  at  the  end.  Again 
the  contour,  is  given  in  Figure  1  and  K+  is  the  plus  factor  of 

K(a)  in  (4)  discussed  in  Appendix  C.  The  integral  may  be  evaluated  ex¬ 
actly  by  finding  the  residues  of  the  poles  of  [K(u)]_1  enclosed  by  Fj 
and  (52)  may  then  be  approximated  by, 


{I  (z)}.  =  V  U  (z.)  I  T  (n)  e 

s°°  int  0  s  0  u ,  On 

n=  1 


(53) 


where  we  have  replaced  the  exact  incident  current,  I^ZgjO)  with  the 
approximate  quantity,  U^z^)  from  (20).  and  have  been  de¬ 

fined  in  (44)  and  (45),  respectively.  Again,  as  in  the  receiving  formula¬ 
tion,  we  may  approximate  the  summation  in  (53)  with  the  expression  in  (47 j. 

It  should  be  noted  that  there  would  also  be  internal  wall  currents 
on  the  semi-infinite  transmitting  antenna  which  would  not  come  from  pene¬ 
tration  at  the  cylinder  end  but  rather  would  be  excited  directly  by  the 
source.  For  a  delta  function  voltage  source,  this  internal  current  can 
be  shown  to  possess  a  logarithmic  singularity  at  the  feed-point  similar 
to  the  logarithmic  singularity  of  the  external  current  at  the  input.  For 
a  more  realistic  excitation,  such  as  a  finite  gap,  however,  the  internal 
current  would  be  well-behaved  everywhere  and  would  be  directly  related  to 

a  capacitive  susceptance  component  (assuming  there  to  be  negligible  radiation 

— y  2 

from  the  open  end  of  the  cylinder  which  in  turn  implies,  e  <<  1) 

of  the  overall  input  admittance.  And  since  we  have  not  addressed  ourselves 
to  the  task  of  specifically  defining  a  "realistic"  input  susceptance,  the 
internal  current  in  the  vicinity  of  the  voltage  source  and  its  effect  on 
the  input  admittance  will  not  be  pursued  any  further  in  this  report. 
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6.3  Internal  current  on  a  f in i te-l.ngth  receiving  antenna 
From  the  external  receiving  current  expression  in  (35),  we  may  write 
the  internal  penetrating  current  near  z  -  -h  using  the  transmission 
characterization  in  (43)  as. 


Iint(0i;z)  =  l  lV(T-0^+h)  Tnn  (^  +  ( n-0 . )  Tn»} 


Y0n(1,+z) 


n=  1 


On 


On 


;  z  -  -h 


(54) 


while  the  penetrating  current  near  the  opposite  end  at  z = +h  is  obtained 
by  replacing  0  with  ir-0^  and  (h+z)  with  (h-z)  in  the  above  ex¬ 
pression.  The  first  term  in  the  {brackets}  above  corresponds  to  TM. 

On 

mode  currents  on  the  internal  walls  of  the  cylinder  due  to  the  primary 

current  term,  E'j'  V(Ti-0.;+h),  while  the  second  term  corresponds  to  TM 

u  l  On 

mode  currents  due  to  the  total  external  current  incident  upon  z  = -h 
arising  from  reflections  eminating  from  the  opposite  end  at  z  = +h. 

These  latter  currents  are  analogous  to  waves  incident  at  an  angle, 

9^  =  7t.  Hence,  the  transmission  coefficient  for  these  incident  currents 
is  evaluated  at  0^  =  n.  The  expression  in  (47)  may  be  used  to  approxi¬ 
mate  the  summations  in  (54),  thereby  reducing  the  computational  efforts 
required  to  find  the  internally  penetrating  current.  Note,  it  is  implicit 
in  this  formulation,  that  there  is  no  internal  interaction  between  the 

”Yoi2h 

ends  of  the  cylinder,  thus  implying  that  e  <<  1. 

6.4  Internal  current  on  a  finite-length  transmitting  antenna 
The  internal  penetrating  current  near  the  ends  of  a  finite  length 
(-h  <  z  <  h)  cylindrical  transmitting  antenna  may  be  obtained  by  applying 
the  transmission  characterization  in  (53)  to  the  respective  incident 
currents  from  (38)  with  the  result. 
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i.  (z  *  +  h)  =  V'  {i;  (h 
int  0  s 


*q)  *  C  (h 


zo)} 


dO 


l 

n=  1 


T0n(7,) 


-y  ( h  ?.  z ) 

On 


(55) 


Analogous  to  the  receiving  case,  the  rirsL  term  in  (55)  corresponds  to 

TM„  mode  currents  on  the  internal  walls  of  the  cylinder  due  to  the 

primary  current,  Vo  Us  ^ '  z  "  Z0 ^  aC  z  =  fh  while  the  second  term 

corresponds  to  TM„  mode  currents  due  to  the  total  external  current 

incident  upon  z =  +h  arising  from  reflections  eminating  from  the  opposite 

end.  Again,  (47)  may  be  used  to  approximate  the  summations  in  (.35)  and 
"Y01(h±z0) 

the  restriction,  e  <<  1,  is  also  implied  in  this  formulation. 


6.5  End  conductance  of  a  finite  length  cylindrical  antenna 
A  quantity  related  to  the  internally  penetrating  current  on  a  cylin¬ 
drical  antenna  is  the  input  conductance  for  a  TM  mode  incident  upon 

On 

one  of  the  ends  of  the  antenna.  Unlike  the  cases  treated  by  Weinstein 

[2,  Chap,  l].  Levin  and  Schwinger  [23],  Jones  [24]  and  others,  our  analysis 

for  the  end  conductance,  discussed  in  Appendix  F,  deals  with  a  TM  mode 

*  On 

under  cut-off,  the  radiation  in  this  case  necessarily  coming  from  tunneling. 
The  end  conductance  in  this  situation  is  relevant  and  very  important  to 
EMC  studies  involving  the  penetration  into  the  end  of  a  cylindrical  en¬ 
closure  [25].  A  detailed  discussion  of  this  quantity  is  left  to  Appendix 

F,  where  the  end  conductance  as  seen  by  an  evanescent  TM„  mode  inside 

Un 

and  near  the  end  of  a  finite  length  cylinder  based  upon  Wiener-Hopf 


analyses  and  the  multiple  reflection  concept  is  derived. 
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Section  7  Numerical  results  for  the  finite  length  cylindrical  antenna 

Due  to  the  restriction  our  theory  places  upon  the  electrical  radius 
of  ka  <  1  ,  the  currents  on  the  internal  wall  of  the  finite  length 
cylindrical  receiving  or  transmitting  antenna  are,  in  general,  very  much 
smaller  in  magnitude  than  the  currents  on  the  external  wall,  except  in 
the  near  vicinities  of  the  ends.  From  Section  6,  it  may  be  ascertained, 
that  the  internal  current  is  significant  only  within  a  distance,  2a  , 
(equal  to  one  cylinder  diameter)  from  either  end.  And  since  we  cannot 
rely  upon  results  from  our  external  current  expressions  so  close  to  the 
ends,  where  the  internal  currents  are  significant,  the  formation  of  total 
current  distributions  from  the  combinations  of  our  receiving  and  trans¬ 
mitting  external  current  distributions  in  (35)  and  (58),  respectively, 
with  the  corresponding  internal  current  distributions  in  (54)  and  (55) 
would  be  of  little  advantage.  Hence,  in  most  cases,  the  external  current 
formulas  in  (35)  and  (38)  will  be  sufl'icien  .  to  describe  the  current 
distribution,  whether  it  be  the  total  or  external  only,  on  finite  length 
receiving  or  transmitting  antennas,  respectively.  On  the  other  hand,  the 
internal  current  distributions  given  in  (54)  and  (55)  for  the  finite 
length  receiving  and  transmitting  cylindrical  antennas  may  be  accurately 
calculated  using  the  approximate  formula  in  (47)  at  practically  ny  point 
on  the  antenna. 

7.1  Current  distribution  on  a  receiving  antenna 

In  order  to  examine  differences  in  our  receiving  theory  resulting 
from  the  use  of  either  approximate  form  of  tJ(9  ;z)  in  (9)  or  (13)  we 
have  included  Figure  17  which  shows  the  magnitude  of  the  induced  current 
at  the  center  of  a  receiving  antenna  where  fl(h)  *  10  ,  illuminated  by 


(A/duj) 


Lectrical  half-length,  kh,  and  normalized  to  the  incident  electric  field 


a  normally  incident  (0  =  tt/2)  plane  w  ,ve  poLarized  parallel  to  the 

antenna  as  a  function  of  the  electrical  length,  using  both  (9)  and  (13) 
in  the  finite  length  receiving  antenna  current  expression  in  (35).  As 
expected,  the  agreement  between  both  results  is  very  good  except  near 
resonances  and  anti-resonances.  And  comparisons  with  existing  analytical 
and  numerical  results  for  cylindrical  antennas  in  which  the  condition, 
(1(h)  >>  | £n(Vg) |  in  (8),  is  satisfied,  have  indicated  that  our  theory 
yields  slightly  better  results  in  almost  every  case  when  (13)  is  used 
for  U(0^;  ?) .  For  these  reasons,  in  what  follows  we  shall  present  only 
results  obtained  from  the  use  of  (13)  for  U(0^;z)  in  the  receiving 
and  transmitting  expressions  in  (35)  and  (38),  respectively. 

The  current  distributions  on  a  half-wave,  kh  =  tt / 2  ,  receiving 
antenna  where  (1(h)  =  2  &n(2h/a)  =  10  for  the  incident  angles, 

0^  =  tt/36,  tt/6,  tt/3  and  tt/2  as  calculated  from  (35)  are  shown  in 
Figure  18.  For  comparison,  first  order  results  from  the  King-Middleton 
theory  [11,  Chap.  IV,  Sec.  7]  and  results  from  King's  three  term  theory 
[5]  for  the  normal  incidence  case,  0^  =  tt/2  ,  are  also  shown.  The 
agreement  between  the  latter  King  theory  and  ours  in  this  particular  case 
is  excellent.  And  the  overall  agreement  between  all  theories  is  quite 
acceptable.  We  note  that  in  spite  of  the  condition  in  (8)  which  requires 
(1(h)  =  2  £n(2h/a)  »  |£n(vg)j  ,  the  current  distribution  predicted  by  our 
formulas  in  the  near-grazing  situation,  0  =  n/36  ,  is  at  a  physically 

anticipated  small  level.  This  is  further  exemplified  in  Figure  19,  where 
for  the  same  antenna  as  in  Figure  18  the  currents  at  z  =  0  ,  h/3  and 
2h/3  are  illustrated  as  a  function  of  the  incident  angle,  0^  .  And  we 
note  the  near  sinusoidal  variation  of  the  current  with  respect  to  the 
incident  angle,  0^  ,  as  would  be  expected  for  a  thin  half-wave  dipole. 
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Figure  18. 


Current  distribution  on  a  thin  half-wave  receiving  antenna 


normalized  to  the  incident  electric  field  and  the  wavelength, 
i.e.,  IR(0.;z)/XE^  . 


from  eq.  (35)  with  (13) 
used  for  U(6  ;z) 


.  0^  =  tt/2  Three  term  theory  of  King  [5] 

.  9,  =  n/2  First  order  King-Middleton  theory 

1  [11,  Sec.  IV. 7] 


Figure  19.  Current  at  specific  positions  on  a  thin  half-wave  receiving 

antenna  as  a  function  of  the  incident  angle,  0^,  calculated 

from  eq.  (35)  with  (13)  used  for  IKG^jz),  and  normalized 

to  the  incident  electric  field  and  the  wavelength,  i.e., 

^  ,>„i 


The  currents  on  the  internal  walls  of  i  lie  same  receiving  antenna 

(kh  =  n/2  ,  (2(h)  =  10  )  at  the  end,  z  -  -h  ,  and  slightly  away  from  the 

end,  z  =  -h  +  2a  ,  as  calculated  from  (54)  with  (47)  are  shown  in 


Figure  20  as  a  function  of  the  incident  angle,  6^  .  The  internal  current 

at  the  end,  z  =  -h  ,  represents  an  infinite  summation  of  all  the  TM 

mode  currents  at  this  point  and  is  equal  in  theory  to  the  negative  of  the 

external  current  at  this  end.  While  the  internal  current  at  z  =  -h  +  2a 

is  predominantly  associated  with  the  TM^  circular  waveguide  mode,  all 

the  higher  order  modes  being  much  more  attenuated  at  this  point.  Thus 

beyond  z  =  -h  +  2a  ,  the  internal  current  will  decay  essentially  as 
-Y01(h+z) 


The  current  distribution  on  an  electrically  thick  (ka  =  1.0) 
receiving  antenna  three  wavelengths  in  length  as  calculated  from  (35)  is 
shown  in  Figure  21  for  the  incident  angles,  0^  =  tj/36  ,  tt/6  ,  tt/3  ,  and 
tt/2  .  Note  that  this  distribution  corresponds  only  to  the  external 
azimuthally  uniform  z-directed  current  on  the  cylinder.  Note  also  that 
since  (2(h)  =  2  Jln(2h/a)  =  5.87  and  j  ?,n(Vg)  j  =  2.64,  0.93,  2.24,  and 
3.63  for  the  respective  angles  considered,  the  condition  that 
(2(h)  >>  |fl,n(Vg)|  as  originally  required  in  the  analytical  development, 
no  longer  holds.  However,  the  correspondence  with  the  data  from  Wu,  et  al. 
[35]  based  upon  the  integral  equation  and  product  integration  formulation 
of  Kao  [ 1 8 ]  for  the  azimuthally  uniform  z-directed  current  also  shown  in 
Figure  21  for  the  same  antenna  with  a  normally  incident  plane  wave  is 
surprisingly  good.  Again  we  bring  attention  to  the  relatively  small  level 
of  current  on  the  antenna  predicted  by  our  theory  at  near-grazing 
incidence,  0^  ■  ir/36.  The  behavior  of  the  current  at  z  *  0  ,  h/3  and 
2h/3  with  respect  to  the  incident  angle,  9^  ,  is  shown  in  Figure  21  and 


INTERNAL  CURRENT  (RECEIVING  ANTENNA) 


RECEIVING  ANTENNA  CURRENT 


kh  =37r 


1  5“ 
u;  4  - 

t  3-^ 

i  2#\  / 


A  /  \ 

; '  i\ 
ft  i  \  /  V  ' 

A  A  A 


j  v 


I  Nl/ 


r  v/ 


-l.O  -0.8  -0.6  -0.4 -0.2  0  0.Z  0.4  0.6  0.8  1.0 

°[  — 1 — H 71  F 

-90°  -  - 

•i8o,k^\/AA; 

-97n«_  *-  y  ^  ^  ^ 


uj 

CO 

<  -270° 
x 

-360° 


Figure  21.  Current  distribution  on  a  three  wavelengths  long,  electrically 

thick,  receiving  antenna  normalized  to  the  incident  electric 

R  i 

field  and  the  wavelength,  i.e.,  1  (S^jzJ/XEq  . 


0J  *  tt/2 
»  it/6 
»  n/3 
-  n/36 
0A  *  tt/2 


from  eq.  (35)  with  (13) 
used  for  U (©^ ; z) 


Numerically  evaluated  data 
from  Wu,  et  al.  [35] 


Figure  22.  Current  at  specific  positions  on  an  electrically  thick 

receiving  antenna  as  a  function  of  the  incident  angle,  0^, 
calculated  from  eq.  (35)  with  (13)  used  for  U(0^;z),  and 
normalized  to  the  incident  electric  field  and  the  wavelength, 
i.e.,  IR(0 . ;z)/AEq  . 
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Figure  23.  Internal  current  at  the  end  and  slightly  within  the  end 
on  an  electrically  thick  receiving  antenna  as  a  function 
of  the  incident  angle,  i)^,  ol  the  uniform  plane  wave, 
calculated  from  eq.  (35)  with  (13)  used  for  U(0^;z),  and 
normalized  to  the  incident  electric  field  and  the  wavelength 
t  e  -  • 


is  seen  to  exhibit  the  physically  expected  result  of  zero  current  ac 
grazing  incidence,  6^=0  and  if  .  The  currents  on  the  internal  walls 
of  this  receiving  antenna  (kh  =  in  ,  ka  *  1)  at  the  end,  z  =  -h  ,  and 
slightly  away  from  the  end,  z  =  -h  +  2a  ,  as  calculated  from  (54)  with 
(47)  are  shown  in  Figure  23  as  a  function  of  the  incident  angle,  0^  . 
Comments  similar  to  the  ones  given  for  the  internal  currents  illustrated 
in  Figure  20  are  also  applicable  to  this  much  thicker  and  longer  antenna. 

7.2  Current  distribution  on  a  transmitting  antenna 

As  discussed  at  the  beginning  of  this  section,  the  total  (internal 
+  external)  current  distribution  on  those  cylindrical  antennas  (both 
transmitting  and  receiving)  for  which  out  theory  is  applicable  is  for  all 
practical  purposes  given  by  the  external  current  distribution  alone, 
except  in  the  near  vicinity  of  the  ends.  An  additional  exception  to  this, 
which  is  particular  to  the  transmitting  antenna,  is  the  region  very  close 
to  the  source  where  internal  currents  are  directly  excited  by  the  source, 
itself.  A  brief  disucssion  of  this  localized  internal  current  has 
already  been  given  in  Section  6.4,  where  it  was  deemed  inappropriate  to 
pursue  an  in  depth  study  of  this  current,  which  is  of  secondary  importance. 

The  current  distribution  on  a  center-driven  half-wave  antenna  where 
fl(h)  *  2  £n(2h/a)  *  10  as  calculated  from  (38)  is  shown  in  Figure  24 
along  with  corresponding  data  from  the  three-term  theory  of  King  [5]  and 
the  approximate  second  order  iteration  procedure  of  King  and  Middleton 
[11,  Chap.  II,  Sec.  22].  The  agreement  between  our  results  and  the  latter 
theory  with  regard  to  the  real  component  of  the  current  is  excellent. 

And  although  the  agreement  between  the  imaginary  components  is  acceptable, 
the  discrepancy  here  was  not  totally  unexpected  since  in  the  process  of 


Figure  24.  Current  distribution  on  a  Lliin,  center  fed,  half  wave 
transmitting  antenna. 

-  from  eq.  (38)  with  (13)  used  for  U(ir  ;z). 

-  -  approximate  second  order  King-Middleton  theory 

[11,  Chap .  II.  Sec.  22] 


King  thrce-tei  i  theory  [5] 
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achieving  an  accurate  value  ^pr  the  real  component  of  the  primary  current 
discussed  in  Section  3.2,  a  less  accurate  "physically  acceptable"  value 
of  the  imaginary  current  near  the  source  resulted.  The  correspondence 
between  the  three-term  theory  of  King  (which  may  be  judged  to  be  less 
accurate  [5]  than  the  King-Middleton  results)  and  our  theory  is  also 
quite  acceptable. 

Perhaps  more  important  than  the  transmitting  current  distribution, 
is  the  input  admittance  to  the  antenna.  Therefore,  in  Figures  25  and  26 
we  show  the  input  conductance  and  susceptance,  respectively,  as  calcu¬ 
lated  from  (41)  for  a  center-driven  cylindrical  antenna  where 
ft(h)  *  2  Hn(2h/a)  =  10  ,  as  a  function  of  the  electrical  length,  kh  . 
Corresponding  admittance  data  from  the  three-term  theory  of  King  [5]  and 
second  order  results  from  the  iterative  method  of  King  and  Middleton 
[11,  Chap.  II,  Sec.  30]  are  also  shown  in  these  figures.  The  agreement 
between  the  conductances  predicted  by  all  three  theories  in  Figure  25  is 
seen  to  be  very  good.  The  agreement  between  the  input  susceptances  is 
also  very  good  for  the  smaller  values  of  kh  where  the  "realistic" 
imaginary  component  of  the  primary  current  is  small  compared  to 
Imaginary  current  arising  from  the  multiple  reflections  from  the  ends. 

At  the  larger  values  of  kh  ,  where  ka  is  proportionally  larger,  we 
find  larger  discrepancies  between  our  results  and  the  King  three  term 
and  King-Middleton  results.  Here  the  imaginary  component  of  the  primary 
current  significantly  affects  the  overall  input  susceptance.  And  since 
our  approximate  expression  for  the  primary  transmitting  current  in  (20) 
is  not  expected  to  accurately  estimate  the  "realistic"  value  for  the 
imaginary  input  current,  this  discrepancy  will  also  appear  in  the  finite 
length  antenna  susceptance  calculated  from  (41)  in  which  (20)  is  used. 
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However,  in  practical  situations  the  input  susceptance  may  be  eliminated 
by  appropriate  matching  leaving  the  input  conductance  essentially 
unchanged  and  the  most  important  quantity  of  consideration. 

To  provide  further  comparison  of  our  theory  with  existing  approaches 
we  offer  Figure  27  which  shows  the  input  conductance  to  a  center-driven 
cylindrical  antenna  where  the  ratio  of  antenna  half-length  to  radius  is 
h/a  =  100  as  calculated  by  (41)  and  the  corresponding  numerically 
evaluated  (via  the  moment  method)  results  of  Harrington  and  Mautz  [26]. 

The  agreement  between  our  results  and  the  accurate  numerically-determined 
data  is  excellent.  Further  evidence  to  substantiate  our  theory  is  given 
in  Figure  28,  which  is  the  same  as  the  previous  figure  except  the  driving 
point  is  now  located  at  z  =  ±h/2  .  Excellent  correspondence  with  the 
numerically  determined  data  of  Harrington  and  Mautz  [26]  is  once  more 
attained. 

We  extend  our  considerations  to  much  thicker  antennas  with  Figure 
29,  which  shows  the  input  conductance  of  cylindrical  antennas  as  calcu¬ 
lated  from  (41)  for  the  radii  normalized  to  wavelength,  a/X  =  0.0159, 
0.078,  and  0.164  (ka  =  0.1,  0.49  and  1.03  ,  respectively)  as  a  function 
of  the  normalized  half-length,  h/A  ,  between  0.1  and  0.5  .  And 
although  these  antennas  are  out  of  the  applicable  range  of  our  theory  due 
to  the  basic  condition  in  (8),  we  find  behavior  still  consistent  with  the 
numerically-obtained  results  of  Chang  [27]  and  [28]  (one-sided  delta 
function  excitation  data),  and  the  experimental  results  of  Hartig  [29]. 

This  further  enhances  the  feeling  that  the  derived  result  actually  has  a 
much  wider  application  than  had  been  assumed  analytically.  And  we  note  the 
very  good  agreement  between  our  theory  and  the  others  for  the  larger  values 
of  h/X  ,  where  the  ratio  of  h/a  is  also  larger. 
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Figure  27.  Input  conductance  of  a  center  fed  cylindrical  antenna  as 
a  function  of  the  normalized  length,  2i  \  . 


However,  in  practical  situations  the  input  susceptance  may  be  eliminated 
by  appropriate  matching  leaving  the  input  conductance  essentially 
unchanged  and  the  most  important  quantity  of  consideration. 

To  provide  further  comparison  of  our  theory  with  existing  approaches 
we  offer  Figure  27  which  shows  the  input  conductance  to  a  center-driven 
cylindrical  antenna  where  the  ratio  of  antenna  half-length  to  radius  is 
h/a  =  100  as  calculated  by  (41)  and  the  corresponding  numerically 
evaluated  (via  the  moment  method)  results  of  Harrington  and  Mautz  [26]. 

The  agreement  between  our  results  and  the  accurate  numerically-determined 
data  is  excellent.  Further  evidence  to  substantiate  our  theory  is  given 
in  Figure  28,  which  is  the  same  as  the  previous  figure  except  the  driving 
point  is  now  located  at  z  =  ±h/2  .  Excellent  correspondence  with  the 
numerically  determined  data  of  Harrington  and  Mautz  [26]  is  once  more 
attained. 

We  extend  our  considerations  to  much  thicker  antennas  with  Figure 
29,  which  shows  the  input  conductance  of  cylindrical  antennas  as  calcu¬ 
lated  from  (41)  for  the  radii  normalized  to  wavelength,  a/A  =  0.0159, 
0.078,  and  0.164  (ka  =  0.1,  0.49  and  1.03  ,  respectively)  as  a  function 
of  the  normalized  half-length,  h/A  ,  between  0  1  and  0.5  .  And 
although  these  antennas  are  out  of  the  applicable  range  of  our  theory  due 
to  the  basic  condition  in  (8),  we  find  behavior  still  consistent  with  the 
numerically-obtained  results  of  Chang  [27]  and  [28]  (one-sided  delta 
function  excitation  data),  and  the  experimental  results  of  Hartig  [29]. 

This  further  enhances  the  feeling  that  the  derived  result  actually  has  a 
much  wider  application  than  had  been  assumed  analytically.  And  we  note  the 
very  good  agreement  between  our  theory  and  the  others  for  the  larger  values 
of  h/A  ,  where  the  ratio  of  h/a  is  also  larger. 


Figure  28.  Input  conductance  of  a  off-center  fed  cylindrical  antenna 
as  a  function  of  the  normalized  length,  2h/A  . 
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The  application  of  our  theory  to  the  thicker  cylindrical  antennas  is 
more  successful  when  the  antenna  length  is  significantly  larger  than  the 
antenna  radius.  Hence,  the  input  admittance  to  a  center-driven  cylin¬ 
drical  antenna  having  a  half-length  to  radius  ratio  of  h/a  =  10  as 
calculated  from  (41)  is  shown  in  Figure  30,  as  a  function  of  the  nor¬ 
malized  length,  2h/A  .  Comparison  of  this  data  to  the  numerically 
obtained  results  of  Harrington  and  Mautz  [26]  also  included  in  the 
figure  is  very  good,  even  near  2h/A  =  2.0  where  the  electrical  radius 
approaches  ka  =  1.256.  This  value  of  electrical  radius  is  slightly 
beyond  the  range  of  our  formulas,  however,  because  the  attenuation  rate 
of  the  secondary  currents  on  an  electrically  thick  cylinder  is  so  pro¬ 
nounced,  as  may  be  observed  in  Figure  14,  the  input  conductance  for  the 
finite  length  electrically  thick  antenna  is  predominantly  determined  by 
Us(0)  in  (20),  which  obviously  still  predicts  the  input  conductance  of 
an  infinite  cylindrical  antenna  to  a  sufficient  degree  of  accuracy. 

Further,  we  offer  Figure  31,  which  is  similar  to  the  previous 
figure  except  the  feed  point  is  taken  to  be  at  z^  =  ±h/2  .  Again  the 
agreement  with  the  numerical  moment  method  results  of  Harrington  and 
Mautz  [26]  is  quite  good. 
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Figure  30.  Input  conductance  of  a  center  fed  cylindrical  antenna  as 
a  function  of  the  nonnnlized  length,  2h/X  . 
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Section  8  Electrically  short  cylindrical  antennas 

8.1  Short  receiving  antenna 

The  current  distribution  on  an  electrically  short  (kh  =  0.4)  thin 
((2(h)  »  2  £n(2h/a)  =  10)  cylindrical  receiving  antenna  calculated  from 
(35)  with  (13)  is  shown  in  Figure  32  for  the  incident  angles  of  the 
uniform  plane  wave  of  8^  =  tt/36  ,  tt/6  ,ir/3  ,  and  tt/2  .  For  comparison 
the  current  distribution  for  normal  incidence  (0^  =  tt/2)  predicted  by 
the  short  antenna  theory  of  King  [11,  Chap.  IV,  Sec.  7]  is  also  shown. 

The  agreement  between  the  two  normal  incidence  magnitude  distributions 
is  quite  acceptable,  especially  in  light  of  the  fact  that  the  basic  con¬ 
dition  in  (8)  is  only  moderately  satisfied,  i.e.,  (2(h)  =  10  while 

| £n(2kh) j  =  0.92  .  On  the  other  hand,  the  phase  of  the  current  predicted 
by  our  theory  is  less  than  90°  along  most  of  the  antenna  as  opposed  to  the 
phase  predicted  by  King  which  is  always  slightly  more  than  90°.  However, 
since  the  phase  of  the  receiving  current  is  quite  inconsequential  compared 
to  the  magnitude  in  most  applications,  this  slight  inconsistency  is  not 
seen  as  a  serious  drawback  to  our  receiving  theory. 

The  behavior  of  the  current  at  the  positions  z  =  0  ,  h/3  ,  and 
2h/3  on  the  same  electrically  short  (kh  =  0.4  and  (2(h)  *  10)  receiving 
antenna  just  discussed  is  shown  in  Figure  33  as  a  function  of  the  incident 
angle,  0^  .  The  sinusoidal  variation  of  the  current  at  each  position, 
with  respect  to  the  angle,  6^  ,  is  readily  apparent  in  this  figure  and 
is  consistent  with  the  expected  behavior. 

8.2  Short  transmitting  antenna 

Figure  34  illustrates  the  current  distribution  on  an  electrically 
short  (kh  ■  0.4)  thin  (J2(h)  -  2  £.n(2h/a)  ■  10)  center  driven  cylinder 
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Figure  33.  Current  at  specific  positions  on  an  electrically  short 

receiving  antenna  as  a  function  of  the  incident  angle,  0^, 
calculated  from  (35)  with  (13)  usefl  for  U(0.;z)  and  normalized 
to  the  incident  electric  field  and' the  wavelength,  i.e., 

IR(9 .;z)/XEj  . 
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Figure  34.  Current  distribution  on  a  center-fed  electrically  short 
cylindrical  transmitting  antenna. 
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calculated  from  (38).  Also  shown  in  this  figure  is  the  current  distri¬ 
bution  on  the  same  antenna  as  predicted  by  the  short  antenna  theory  of 
King  [11,  Chap.  II,  Sec.  31].  The  agreement  between  our  result  and  that 
of  King  is,  in  general,  fairly  good,  especially  since  we  have  again  only 
moderately  satisfied  the  basic  restriction  in  (8).  Of  particular 
interest  here  is  the  inability  of  our  formulation  to  predict  an  accurate 
real  component  of  current  at  the  source.  This  situation  results  from 

the  fact  that  our  transmitting  current  formulation  in  (38)  is  good  only 
-2 

to  (C^)  for  short  antennas,  while  the  input  conductance  for  these 
antennas  is  given  by  higher  order  terms.  We  may  easily  recover  the  input 
conductance,  however,  by  turning  to  the  input  conductance  formulation  of 
Chang  and  Rispin  in  [30]  which  is  based  upon  the  effective  aperture  of  the 
antenna  and  requires  a  knowledge  of  only  the  magnitude  of  the  receiving 
current  on  the  same  antenna  illuminated  by  a  normally  incident  uniform 
plane  wave.  Their  approximate  result  for  the  electrically  short  receiving 
antenna  may  be  written  as 


G 


j  I R(tt/ 2 


(56) 


where  |IR(tt/2  ,  z)/Eq|  is  the  magnitude  of  the  receiving  current  nor¬ 
malized  to  the  incident  field,  at  the  feedpoint,  Zq  ,  due  to  a  normally 
incident  (0^  =  tt/ 2)  plane  wave.  As  seen  in  the  previous  sub-section, 
the  magnitude  of  our  receiving  current  is  of  sufficient  accuracy  for  these 
short  electrical  lengths  to  permit  the  above  calculation.  Figure  35  shows 
the  input  conductance  predicted  by  (56)  using  our  receiving  current 
formula  in  (35)  along  with  comparable  data  from  the  short  antenna  theory 
of  King  [11,  Chap.  II,  Sec.  31]  for  electrically  short  antennas  charac¬ 
terized  by  Q(h)  =  2  2,n(2h/a)  *  10  as  a  function  of  the  electrical 
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Figure  35.  Input  conductance  of  a  center-fed  electrically  short 
cylindrical  antenna. 
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length,  between  kh  =  0.1  and  1.0.  Obviously  the  two  theories  appear 
to  be  in  relatively  good  agreement.  Thus,  even  though  our  transmitting 
current  formula  in  (38)  failed  to  yield  an  acceptable  input  conductance 
for  the  electrically  short  cylindrical  antenna,  the  alternate  procedure 
described  above  that  relies  upon  our  compl ementarv  receiving  current 
formula  in  (35)  has  yielded  successful  results. 
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9.  Concluding  remarks 

Through  a  re-examination  of  the  conditions  necessary  to  obtain  simple 
asymptotic  solutions  based  upon  a  Wiener-Hnpf  analysis  for  the  reflected 
current  distribution  on  semi-infinite  tubular  cylindrical  antennas,  we  have 
broken  away  from  the  traditional  "thin  wire"  assumptions,  ka  <<  1  and 
kh  >  1  ,  and  extended  these  asymptotic.  Wiener-Hopf  solutions  to  include 
antennas  as  thick  as  ka  =  1  and  as  short  as  kh  =  0.4  .  These  extended 
solutions  have  been  derived  subject  primarily  to  the  satisfaction  of  the 
condition,  11(h)  =  2  £n(2h/a)  >>  [  in[2khsin~  (Ch/2)  ]  |  .  Although  for  the 

thicker  antennas,  we  have  observed  a  significant  relaxation  of  tl  >s  condi¬ 
tion.  Our  simple  receiving  current  formula  is  applicable  for  all  angles  of 
incidence  of  a  uniform  plane  wave,  i.e.,  0  <  0.  <  tt  ,  and  we  have  con¬ 
structed  our  simple  transmitting  formula  in  such  a  wav  so  that  an  accurate 
input  conductance  is  obtained  over  our  entire  range  of  interest.  Numerical 
comparisons  of  our  theory  with  the  work  of  many  other  authors  have  shown 
very  good  overall  agreement.  Hence  our  theory,  which  requires  primarily 
only  arithmetic  calculations,  permits  easy  and  inexpensive  calculations  of 
the  current  distribution  on  cylindrical  transmitting  and  receiving  antennas, 
thus  providing  an  excellent  basis  for  statistical  and  transient  studies. 

As  a  consequence,  the  transient  aspects  of  a  cylindrical  antenna  based  upon 
our  present  theory  shall  be  covered  in  a  forthcoming  report. 

Although  only  the  unloaded  transmitting  and  receiving  antenna  has  been 

analyzed  in  this  report,  the  extension  to  the  loaded  situation  in  both  cases 

is  straightforward.  The  voltage  developed  across  a  load  admittance.  Y(  , 

located  at  z  =  z^  along  a  cylindrical  transmitting  or  receiving  antenna 

T  X  R  R 

would  be  given  by  V?  =  -I  ‘(zQ  ;zp) /Y  or  V;i  =  -I  ?  (0  ;  )  /  Y  ^  ,  respectively, 

T  R 

where  I^Uq.z^)  and  1^(6. ;z^)  are  the  transmitting  and  receiving  currents. 
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respectively,  at  z  =  z^  in  the  loaded  situation.  The  potential  difference 

across  z  =  z^  in  either  the  transmitting  or  receiving  situations  acts  as 

a  voltage  generator  for  which 'our  transmitting  antenna  theory  can  be  used  to 

obtain  the  current  distribution.  Adding  this  current  distribution  to  the 

unloaded  current  distribution  yields  the  total  (external)  current  distribu- 

T 

tion  which  may  in  turn  be  solved  at  z  =  Zg  to  obtain  I^Cz^z^)  or 
R 

l^(0^;z^)  .  Subsequently  the  loaded  receiving  and  transmitting  current 
distributions  may  be  written  as  [11,  Sec.  IV. 7], 


and. 


.  lR«yv  ‘byo 

«(  C2)  ‘  CZ)  ’  Y1+Yln(Ze)  V0 


.t,  ,  tt,2  . 

VVz)  "  1  0’Z)  "  Y  +Y  (2  )  Vq 


(57) 


(58) 


respectively.  The  unloaded  receiving  current,  I  ((.y;z)  is  given  in  (35) 

T 

and  the  unloaded  transmitting  current,  1  (z^;z)  in  (38).  The  term, 

T 

I  (z^;z)/Vq,  corresponds  to  the  transmitting  antenna  current  distribution 
due  to  a  voltage  source  of  unit  strengtli  at  z  =  zf  .  And  Y^(Zg)  is  the 
input  admittance  to  the  antenna  at  ::  =  z{  which  is  approximately  given  by 
(41).  Although  (41)  has  been  shown  to  give  an  excellent  result  for  the 
input  conductance,  i.e.,  Re{Y^n)  ,  the  input  susceptance,  -ImfY^}  from 
(41)  is  subject  to  error.  Thus  in  some  cases  a  more  accurate  expression  for 
the  input  admittance  Y^n(zj^)  may  be  desired.  The  error  in  our  admittance 
formula  in  (41)  comes  mainly  from  the  inaccurate  prediction  of  the  input 
susceptance  by  Ug(0)  for  an  infinitely  long  antenna.  For  the  smaller  values 
of  ka  (ka  <  0.1)  the  imaginary  part  of  the  term,  -U  (0)  in  (41),  may  be 
replaced  by  the  susceptance  formula  of  Fante  [21],  which  has  been  shown  to  be 
a  good  approximation  for  the  input  susceptance  by  Miller  [32]  and  may  be 
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written  as. 


£n(kd) 


where  6  is  the  physical  width  of  the  load  region  and  is  assumed  to  be  very 


small.  As  noted  by  Miller  [32],  the  input  conductance  to  an  infinite 
cylinder  is  relatively  insensitive  to  variations  in  the  gap  width,  6  ,  when 
6  is  small  compared  to  the  cylinder  radius,  a  .  Hence,  Re{Us(0)}  may 
be  retained  in  (41)  to  determine  the  infinite  cylinder  input  conductance. 
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The  current  distributions  on  antennas  with  multiple-loads  can  be  found  by 
generalization  of  the  single-load  approach  outlined  here. 

Another  quantity  of  interest  for  which  our  approximate  cylindrical 
antenna  theory  proves  useful  is  the  far  field  radiation  pattern  of  a 
cylindrical  transmitting  antenna.  It  is  well  known  that  the  far  field 
radiation  pattern  of  a  transmitting  antenna  bears  the  same  angular  depen¬ 
dencies  as  the  received  current  at  the  same  feed  point  when  the  antenna  is 
used  as  a  receiving  element.  A  derivation  of  this  principle  as  applied  to 
the  cylindrical  antenna  is  contained  in  Appendix  H.  From  Appendix  H,  we 
have  the  far  field  radiation  from  a  finite  length,  -h  <  z  <  h  ,  cylindrical 
transmitting  antenna  with  a  delta  function  voltage  source  of  strength, 
volts,  at  z  =  Zq  given  in  terms  of. 


r  .  .  +ikn  elkr  ..  1  ;z0> 

E„(r,e.«  -  — —  v0 - - - 


kr  >>  1 


(r,Q,4i)  refers  to  a  spherical  co-ordinate  system  coincident  with  the 

cylindrical  co-ordinate  system  implied  in  Figure  16  ,  i.e.,  the  positive 

R  i 

z-axis  corresponds  to  0=0.  And  I  (it  -  0  jz^/E*  is  the  received  current 
at  z  •  Zq  on  the  same  antenna  when  illuminated  by  a  uniform  olane  wave  of 
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unit  amplitude  incident  at  an  angle,  r-0  ,  with  respect  to  positive 

z-axis.  Hence  the  far  field  radiation  from  a  transmitting  antenna  ma  be 
approximately  determined  from  (60)  using  our  approximate  formula  for  the 
receiving  current  distribution  in  (35)  for  the  same  antenna  parameters. 

For  example,  for  a  center-fed,  Zq  =  0  ,  transmitting  antenna  where 
0(h)  =  2  £n(2h/a)  =  10  and  kh  =  n/2  (the  approximate  current  distribution 
of  which  is  shown  in  Figure  24),  the  far  field  radiation  is  shown  in  the 
z  =  0  curve  (apart  from  the  constant  exp(ikr)  )  of  Figure  19  with 

0^  replaced  by  tt  —  0  Similarly  all  of  the  receiving  antenna  current 
distributions  shown  in  graphical  form  in  this  report  may  be  related  to  the 
far  field  of  the  same  antenna  when  it  is  used  as  a  transmitting  element . 

The  basic  approach  utilized  in  this  report  to  develop  a  simple  approxi  - 
mate  cylindrical  antenna  theory  is  currently  being  applied  to  other  types 
of  linear  antennas,  such  as  the  co-axial  and  parallel  thin-wire  antennas. 
Results  of  this  research  should  appear  in  the  not-too-distant  future. 

In  the  same  respect,  since  our  present  theory  deals  only  with  the  azimuthal ly 
uniform  axial  current  on  cylindrical  antennas,  a  further  extension  of  our 
theory  which  also  considers  the  higher-order  variations  of  axially-directed 
currents  as  well  as  the  circumferentially-directed  currents  is  currently 


being  prepared. 
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Appendix  A.  Approximate  solutions  for  the  canonical  integral.  U(0.;z) 


Before  the  "so-called"  canonical  integral  for  U((h;z)  in  (3)  is 
considered,  let  us  first  of  all  begin  with  the  closely  related  auxiliary 
canonical  integral,  W(0^;z)  .given  in  (6)  and  repeated  here, 

,,  ,  ,  -iaz 

W(0 . 

M  i  J 

'0 


-ik  f 

) . ;  z)  =  -  (1  -  cos  0  . )  7i - 771 - — — 

l  *  n  i  j  p  (k  +  a)  (k  cos  0. 


+  a)K(a) 


da 


where , 


and 


0  <  Z  <  °°,  0  <  0 .  <  7T 

-  —  1  — 


,(1) 


K(a)  =  iiU0(^a)H^  '(f,a) 


(Al) 


(A2) 


/.  2  2  .  fT Ti 

=  /k  -  a  =  i/a  -  k 


(A3) 


K+  is  the  "plus"  factor  of  K(a)  which  is  analytic  and  free  of  zeros  in 
the  upper  naif  a-plane,  (see  Appendix  C) ,  and  the  contour,  Fq  ,  is 
shown  in  Figure  1.  We  ^'<hall  find  approximate  solutions  for  W(0^;z)  in 
(Al)  and  a  similar  integral  for  11(0^;^)  .which  is  stated  later,  that 
require  the  satisfaction  of  the  condition. 


fi(z)  «  2  2n(2z/a)  >>  |£n[2kz  s in^ (0^/2) ] | 


( A4) 


Actually,  the  restriction  in  (A4) ,  which  will  be  referred  to  as  the 
"basic  condition"  on  our  analysis,  expresses  the  most  severe  restriction 
encountered  in  the  derivation  of  the  approximate  solutions  for  the 
canonical  integral,  U(6^;z)  ,  and  the  auxiliary  canonical  integral, 
W(0^;z)  .  Many  of  the  approximations  to  follow  require  lesser  restric¬ 
tions  than  the  one  in  (A4)  and,  therefore,  are  automatically  satisfied 
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when  (A4)  is  true. 

Since  K+(-a)  is  analytic  in  the  Lower  half  a-plane,  combining 
the  contributions  from  both  sides  of  the  branch  cut,  which  is  due  to 
the  term  [K(a)]_1  in  the  integrand,  allows  (Al)  to  be  written  as. 


W(6.;z) 


ikz 
e  1 

f  ^ 

[i  i  1 

k  +  ^ 

Z 

-u 

7TH  J 

u  u-iv^ 

IT  e 

0 

Jo 

-  /u(u  -  iZkz) 
z 

lho 


(1)  I  —  /u(u  -  i2kz7 


i 

h(2)  ! 
Ho 

—  /  u(u  -  i2kz) 
z 

du  (A5) 


where  the  change  of  variable,  a  =  -k  -  iu/z  ,  has  been  utilized.  vQ 
is  defined  by. 


and. 


vQ  =  voPi;z^  =  2kz  sin2(0i/2) 


— tt/4  <  arg(  /u(u  ~  i2kz)  )  <  0 


(Afe) 


(A7) 


Due  to  the  exponential  decay  of  the  integrand  in  (A5) ,  the  meaningful 
range  of  integration  is  limited  to  the  values  of  u  of  order  unity  or 
less.  And  since  the  term  K+(a)  depends  only  upon  the  parameters,  ka 
and  aa  ,  (see  Appendix  C  for  a  discussion  of  K+(a)  )  ,  then  under  the 
basic  condition  in  (A4)  (which  implies  a/z  <<  1  )  the  term 
K+(k  +  iu/z)  appearing  in  (A7)  may  be  accurately  represented  by  the 
first  term  of  its  Taylor  series  expansion  about  u  =  0  ,  i.e.,  K+(k)  . 

This  term  may  then  be  removed  from  the  integrand  in  (A5 )  yielding, 

W(0i; z)  =  K+(k)U(9i;z) 


(A8) 
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where  U(0^;z)  is  the  canonical  integral  given  by. 


ikz  r  00 


U(e.;z) 


i-  °°  r 

i _ i _  — 

J0  Lu 


—  /u(u  -  i2kz) 


/ ,  \  — — — p  +  — — — — - - - - 1-  du  (A9) 

—  /u(u  -  i2kz)  H  —  /u(u  -  i2kz) 

U  2  U  2 


and  is  equivalent  to  the  expression  for  U(0.;z)  in  (3)  of  Section  2. 

The  evaluation  of  U(0^;z)  in  (A9)  requires  numerical  integration, 
since  an  exact  analytical  solution  to  the  integral  involved  ;s  not 
known.  Hence,  we  seek  a  simple  approximate  solution  for  U(fK;z) 
which  possesses  a  wide  range  of  applicability.  To  this  end,  we  assume 
the  satisfaction  of  the  basic  condition  on  U(z)  in  (A4) ,  which  allows 
the  Bessel  functions  in  (A9)  to  be  approximated  by  the  leading  terms 
of  their  respective  small  argument  expansions.  The  resulting  approximate 


expression  for  UCO^jz)  can  be  written  as, 


U(0i;z)  z  m(0;z)  -  pCO^z) 


(A10) 


where. 


uOi;z) 


ikz  r  °°  -u  f 


.  O)  ^-U 

Jo  (u  ~ iv0 


)  fi(z)  -  2y  +itr  -  £n[u(u  -  i2kz)  ] 


Q(z)  -  2y  -  in  -  .’.n[u(u  -  i2kz)  ] 


(All) 


and  ft(z)  is  given  in  (A4) .  The  [bracketed]  term  behaves  essentially 
like  a  slowly-varying  function,  both  magnitude  and  phase,  over  the 
integration  range,  0  <  u  <  1  ,  where  the  integrand  has  its  major 


contribution,  except  of  course  very  near  u  =  0  .  Thus,  provided  v„ 


i 


V 
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is  not  very  small,  we  can  approximate  (All)  up  to  order  [S2 ( z)  by. 


.  ikz 

p(0i;z)  =  -V 


" _ 1 _ 

f2(z)  -  2y  -  £n(2kz)  +  i3u/2 


-iv. 


fi(z)  -  2y  -  £n(2kz)  -  in/2 


(A]  2) 


where  is  the  exponential  integral  of  the  first  kind  defined  in 

[16,  Eq.  5.1.1], 

The  situation  is  somewhat  different  for  the  case  when  U .  =0 

l 

because  the  integrand  now  blows  up  at  u  =  0  .  To  evaluate  mis 

integral  properly,  we  must  retain  the  logarithmic  behavior  of  the 

integrand  near  u  =  0  ,  while  approximating  the  term  £n(u-i2kz)  bv 

£n(-i2kz)  .  This  procedure  can  be  shown  to  lead  to  a  solution  for 

_2 

p(0;z)  good  to  order  [f2(z)j  consistent  with  the  )i(0^;z)  solution. 
Proceeding,  we  obtain. 


U(0;z) 


" _ 1__  _ _ 

z)  -  2y  -  £n(2kz)  +  i 3tt/ 2  -  £n(u) 


_  1 _ _ 

ti(z)  -  2y  -  £n(2kz)  -  in/2  -  £n(u) 


du 


(All) 


which  has  been  previously  given  by  Shen  [6]  whose  subsequent  solution 
was  based  upon  the  observation  that  the  integration  of  a  smooth  function 
weighted  by  an  exponential  decay  over  an  infinite  interval  can  be 
approximated  by  the  unweighted  integration  of  the  same  function  over  a 
judiciously  chosen  integration  range.  The  same  result,  however,  can  be 
obtained  in  a  more  straightforward  manner  if  we  replace  u  ^  by 
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— WVl 

e  dw  and  interchange  the  order  of  integration  in  (A13)  which 

1 

yields, 

_ 1_ 

Q(  z)  ~  2y  +  £n(w)  -£n(2kzt)  + i 3r/2 

1 

Q(z)  -  2y  +  £n(w)  -  £n(2kzt)  -  in/2  J  clt  [ 

j 

Now  since  the  term  |£n  t|  is  typicalLv  small  compared  with  jf)(z)| 
for  t  <  1  ,  except  very  close  to  t  =  0  where  the  integration  is 
rather  insignificant,  a  two  term  Taylor  series  expansion  followed  by 
straightforward  integrations  on  t  yields. 


4(0; z) 


le 


dw 

w 


-t 


4(0;z)  * 


ie 


ikz 


dw 

w 


/ _ _ _ i _ r 

-  2y  +  £n (w)_£n(-i2kz)  +  in 


[ 


Y 


S2(z)  -  2y  +  £n  (w  )-£n(-i2kz)  +  it 


fi(z)  -  2y  +  £n  (w)-£n(-i2kz)  -  in 


1  - 


Y 


Q(z)  -  2y  +  £n  (w)-£n(-i2kz)  -  in 

(A15) 


1 


Now  taking  the  [bracketed]  terms  in  (A13)  as  being  the  first  two  terms 
of  a  Taylor  series  expansion  of  a  function  of  the  form  [1+x]  ^  for 
small  x  ,  the  integrals  in  (A15)  can  be  evaluated  to  order  [$2(z)  ]  "  , 
with  the  result  for  p(0;z)  being  given  by. 


U(0;z) 


ieikz  f  Q(z)_  y  -  i ?t  —  ?.n(-i2kz)  T 

H  |fl(z)  -  y  +  in  -  f.n(-i2kz)J 


(Alb) 
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(A12)  and  (A16)  may  be  combined  to  give  a  solution  for  (A10)  good  to 

__*v 

order  (£2(z)  ]  *  which  is  given  by. 


.  ikz  fi(z)-'y-i'n-iin(-i2kz)+e  E  (-iv  ) 


U(0i;2)  =  ln 


lx  O'  ( 


r  (A17) 


z)  -  y  +  in  -  ?.n  (-i  2kz)  +  e  (-iv^) 


Note  that,  although  this  approximate  expression  for  U(0.;z)  was 

2 

derived  under  the  assumptions  that  v^  =  2kz  sin  (0^/2)  cannot  be  too 
small  and  that  |fi(z)  j  >>  |in(  )  |  ,  the  particular  method  of  combining 
y(0;z)  and  y(0^;z)  into  U(9^;z)  in  (A17)  actually  gives  a  smooth 
function  as  z  -*•  0  . 

Although  UCG^jz)  is  expressed  in  fairly  simple  terms  in  ( A 1 7 i . 

_  O 

a  simpler  form  (also  good  up  to  0[Q(z)]  “  )  may  be  obtained  by  express  in 
the  £n  in  terms  of  an  arctangent,  which  in  turn  may  be  approximated 
by  the  leading  term  of  its  Taylor  series  expansion  when  the  basic 
condition  in  (A4)  is  satisfied.  This  procedure  yields  the  alternate 
formula  for  U(d^;z)  given  by. 


U(0i,z) 


( A 1 8 ) 


y  -  £n(-i2kz)  +  e  E^(-iVp) 
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Appendix  B .  Input  conductance  of  an  i nf in  i  t e ly  long  cy  1  ind er 


The  input  conductance  (which  is  associated  with  radiation  from  the 
external  surface)  of  an  infinitely  long  hollow  cylindrical  antenna  may  be 
formulated  by  means  of  considering  the  real  part  of  the  source  current  for 
a  delta  function  voltage  source  of  unit  strength  (i.e.,  1.0  volt).  [22]. 
The  resulting  expression  may  be  written  as. 


G  (ka) 

CO 


Aka 

Tin 


ka 
•  0 


dx 


xv(ka) ' 


x2  J0(X)  +  Y0(x) 


(1)1) 


where  n  is  the  intrinsic  impedance  of  the  surrounding  medium.  is  the 

cylinder  radius,  and  J_  and  Yq  are  the  zero  order  Bessel  and  Neumann 
functions,  respectively.  Although  several  asymptotic  solutions  (which 
assume  ka  <<  1)  for  Gm(ka)  exist,  [3],  [A],  [20]  and  [21],  an  exact 
solution  to  the  integral  in  (Bl)  is  not  known  to  the  authors.  However,  in 
their  Wiener-Hopf  approach  to  thin  dipole  antennas  Shen.  Wu,  and  King  [6] 
constructed  the  surprisingly  accurate  formula  for  Go(ka)  given  by. 


G  (ka)  =  Re  { —  [  £n(  2C  -  iir/2)  -  Cn(2C  +  i3'i/2)]}  (B2) 

°°  n  w  w 

which  agrees  with  the  asymptotic  expression  of  Wu  [3]  to  order  C  '  , 

w 

(C^  =  -Hn(ka)  -  y  ,  y  =  0.577...)  .  However,  the  constructed  expression  for 

Goo(ka)  in  (B2)  is  fairly  accurate  even  when  J C  j  is  not  large.  In 

Figure  Bl,  we  have  illustrated  the  input  conductance  predicted  by  botu 

Shen’s  [6]  curve-fitting  formula  and  the  exact  expression  in  (.Bl)  over  tin 

-A 

range  of  electrical  radii,  10  <  ka  <  1.0  .  This  l igure  reveals  tin- 

remarkable  ability  of  Shen's  approximate  formula  to  predict  an  .nee,:  , 
input  conductance  not  only  in  the  conventional  th in-wire  regions 


■ 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  Of  STANDARDS- 1963-A 


;ure  Bl.  Input  conductance  of  an  infinitely  long  tubular  cylinder. 
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for  cylinders  as  thick  as  ka  ~  0.7  ,  where  C  =  -0.22  !  Furthermore,  as 

w 

ka  -*■  0  ,  the  approximate  formula  given  in  (B2)  becomes  asymptotic  to, 

Goo<ka)  ~  J  (Cw)“1  ;  as  ka  +  0  (B3) 

which  corresponds  to  the  asymptotic  behavior  of  the  results  of  many  other 
authors,  [3],  [4],  [20]  and  [21]. 

In  order  to  further  extend  the  applicable  range  of  the  basic  approxi¬ 
mate  formula  in  (B2)  to  include  values  of  ka  as  high  as  unity,  we  have 
employed  a  curve-fitting  procedure  to  produce  the  modified  formula, 

Gjka)  s  Re  [  An(2Cw  -  g  -  m/2)  -  An(2Cw  -  g  +  i3ir/2)  ] }  (B4) 

where 

g  -  33.88(ka)2  exp  [-  ^  ]  (B5) 

Obviously  the  insertion  of  g  in  (B5)  into  the  approximate  expression  for 
Goo(ka)  in  (B2)  has  little  effect  below  ka  ~  0.5  ,  however,  above  this 
point,  all  the  way  up  to  ka  =  1.0  ,  this  modification  allows  for  a  quite 
accurate  prediction  for  the  input  conductance  as  can  be  observed  in  Figure  5, 
which  is  discussed  in  the  text  of  this  report. 

Consequently,  in  order  to  achieve  this  accurate  value  of  input  conduc¬ 
tance  from  our  primary  transmitting  current  formula  in  (19)  at  z  «  Zq  ,  we 
follow  the  lead  of  Shen,  et  al.  [6]  and  restate  the  approximate  form  of 
U(ir;z)  in  (9)  in  the  modified  form, 

us<z)  e1KZ  Un[fs(z)  -iTT]  -  £n[fg(z)  +iir]}  (B6) 


where. 


fg(z)  -  2Cw  +  Y  +  iir/2  +  in[kz  +  /( kz)2  +  exp(-2y-2g)  ]  (B7) 

As  given  in  (B6) ,  Ug(z)  corresponds  exacLly  to  Shen's  Eq.  6  in  [6]  except 
for  the  additional  term,  -2g  ,  in  the  square  root,  which  enables  (B6)  to 
yield  a  very  accurate  value  for  the  real  component  of  the  input  current  up 
to  ka  =  1.0  .  And  we  note  that  U(tt;z)  in  (9)  and  the  modified  form 
denoted  by  Ug(z)  in  (B6)  differ  insignificantly  for  kz  >  1  . 
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Appendix  C.  The  kernel,  K(g),  and  the  factorized  functions,  Kt  (g) 
and  K  (a) 

Consider  the  Wiener-Hopf  kernel, 

K(a)  -  i*JoUa)Hf(1)aa)  (Cl) 


where 

£  -  /k2-a2  -  i/a^-k2  (C2) 

A  small  argument  approximate  form  of  (Cl)  which  proves  useful  when  dealing 
with  conventionally  thin  (ka<<i)  cylindrical  antennas  is  given 

by. 

\ 

K(a)  =  2Cw  +  irr  -  2  in  )  >  ka,  aa<<l  (C3) 

where 


Cw  -  -  in(ka)  -  y  ,  y  *  0.57721566  (C4) 

In  the  Wiener-Hopf  procedure,  the  kernel  in  (Cl)  is  factorized  into 
the  form. 


K(a)  -  K+(a)K_(a)  (C5) 

in  which  K+(a)  is  analytic  and  free  of  zeroes  in  the  upper  half  a- 
plane,  -Im(k)  <  Ia(a)  <  ®  ,  and  K  (<*)  is  analytic  and  free  of  zeroes 
in  the  lower  half  a-plane,  <  Im(nc)  <  Im(k)  .  We  assume  a  symmetrical 
factorization  which  gives  the  relationship, 

K_(-a)  -  K+(a)  (C6) 
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Thus  the  "minus  factor",  K_(a),  is  known  once  K+(a)  is  found.  Another 
property  imposed  upon  the  factorization  of  K(a)  is  the  asymptotic  be¬ 
haviors  , 


and 


K+(a)  'a  ;  as  |  a  |  -►  °°  in  the  upper  half  a-plane 

(C7) 


K_(a)  -  a  ;  as  |  a  |  ->-  00  in  the  lower  half  a-plane 

(CS) 


For  the  special  case  where  ka  and  aa  are  much  less  than  unity, 
K+(a)  may  be  approximated  by  the  leading  term  of  Hallen's  exact  formula 
for  K+(a) ,  [1, Sec. 38],  given  by. 


K+°  (a)  =  ^2Cw+iir 


ka,  aa<<l 


(C9) 


where  is  given  in  (C4).  The  superscript  "0"  is  used  here  to  signify 

this  as  being  the  small  argument  (i.c.,  thin  wire)  approximate  formula 
for  K+(a)  .  And  we  note  that  (a)  in  (C9 )  may  also  be  identified  as 
the  plus  factor  in  the  approximate  factorization  of  the  small  argument 
approximation  of  K(a)  in  (C3)  given  by: 

K(a)  =  C2S.+1.)  [l  -  ^  >"  (  IT  )] 

;  ka,  aa<<l  (CIO) 

which  is  valid  to  order  (C  )”2  . 

w, 

On  the  other  hand,  the  exact  value  of  K+(a)  may  be  obtained 
numerically  from  the  formula  of  Mittra  and  Lee  f 15, Sec. 5-2. (3)  ]  which 


is  given  by 


K+(a)  -  /U  [jQ(£a)]  +  [H0(1)  (.  a)]+e 


sx(») 


(Cl  1 ) 


(1), 


where  £  is  given  in  (C2).  The  plus  factors  of  J^(£a)  and  (£a) 


are  given  by 


oo 

[J  (£a)]  =  /ykT)'  n  (1  +~  ) 

n=l  '  Y0n 


.  a 

i —  a 


nn 


(C 12) 


and 


[HQ(1)(£a)]  +  =  /hJ1W)  exp  {  -i  y  +  +  q (a)}  (C13) 

respectively.  Here  iyQn  corresponds  to  the  propagation  constant  of  the 


TM.  mode  within  a  circular  waveguide  of  radius  a  and  is  given  by 
On 


nr, 


On 
X(a) 


J\ 


On 


-  k2  ;p  =n  ordered 
n  zero  of  Jq. 


(CIA) 


The  convergence  factor  eA  in  (C 11)  insures  the  correct  asymptotic 
behavior  described  in  (C 7 )  and  is  given  by 


X(a) 


“  t  “[- 


1  +  y  -  In 


/  2l  \ 

_  i  ll 

\  ka  ) 

1 2J 

(C15) 


where  y  m  0.577215 .  And  finally  q(a)  is  given  by 


,(«)-■![  - .1 

*  4  L  X  +  Y02(X)  / (ka)2-x 2  / 


dx 


(C16) 


The  specific  application  of  K+(ot)  in  this  paper  is  to  calculate 
the  reflection  coefficient  R(0.)  ■  (n/2n)K  (k)K  (-k  cos0.)  in  the  range 
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0  £  0^  <  it  ;  thus  we  require  an  accurate  lormula  for  K+(a)  for  the 
corresponding  range  -k  <  a  <  k  .  Such  a  formula  may  be  ob¬ 
tained  by  observing  the  behavior  of  the  ratio  of  the  small  argument 
approximation,  (a)  in  (C8),  with  respect  to  the  exact  form  of  K+(a) 
in  (Cll)  over  the  range  -k  <  a  <  k  .  Since  the  behavior  of  the  small 
argument  formula  could  be  qualitatively  argued  as  being  at  least  partially 
inherent  in  the  exact  formula  for  K+(a)  for  these  values  of  a  ,  the 
ratio,  (a)/K+(a),  would  be  expected  to  be  a  fairly  smooth  function  of 
a  .  This  preassumption  is  verified  in  Figure  Cl  in  which  we  show  the  real 
and  imaginary  parts  of  the  ratio  K®  (a)/K+(a)  for  ka  =  0.01,  0.05, 

0.1,  0.5,  and  1.0,  as  a  function  of  a  between  -k  and  +k  .  Noting 
the  behavior  of  this  ratio  and  its  predominantly  real  nature,  we  may  assume 
the  approximate  relationship. 


_  A[1  4-  B(a/k)]  ;  0  £  a  <  k 

K+°(«)  „  « 

K+  (a)  |_A(1  +  B(a/k)  +  C(a/k)2]  ;  -k  <  a  <  0 


(Cl  7) 


where  A,  B,  and  C  are  coefficients  yet  to  be  determined.  The  above 
linear  relationship  for  0  <  a  <  k  appears  to  describe  the  behavior  of 
the  imaginary  component  of  (a)/K+(a)  for  the  cases  where  ka  <  0. I 
in  Figure  Cl,  while  it  appears  that  higher  order  terms  in  a/k  are 
necessary  above  ka  :  0.1.  However,  due  to  the  relatively  large  linear 
real  component  of  K+  (a)/K+(a)  ,  neglecting  this  higher  order  variation 
in  the  imaginary  component  will  not  produce  a  significant  error.  The 
constant  A  is  found  by  setting  a  -  0  in  (C17)  yielding 
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The  coefficient  B  is  found  by  differentiating  (C17)  with  respect  to  a 
and  setting  a  =  0  ,  the  result  being 
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As  it  stands,  it  appears  as  though  a  numerical  integration  is  necessary 
in  (C19)  to  evaluate  B  .  However,  one  may  recognize  that  the  real  part 
of  the  integral  in  question  is  directly  related  to  the  integral  contained 
in  the  expression  for  the  input  conductance  of  an  infintely  long  cylinder, 
G^(ka)  ,  in  (Bl)  of  Appendix  B.  And  since  a  very  good  approximation  for 
G  (ka),  good  up  to  ka  *  1.0,  has  been  found  and  is  stated  in  (B4),  we 
know  the  real  part  of  B  in  (C19)  to  a  very  good  degree  of  approximation. 
Furthermore,  noting  from  Figure  Cl  that  the  imaginary  component  of 

(a)/K+(a)  for  all  the  values  of  ka  considered  is  relatively  small, 
we  may  conclude  that  the  imaginary  part  of  B  and  the  neglected  higher 
order  imaginary  components  are  not  important.  Specifically  for  the  con¬ 
stant  B  we  may  state  that  , 


2  2 
B.  <<  (l  +  B  ) 
i  r 


(C20) 


where  Bf  and  are  the  real  and  imaginary  parts  of  B  in  (C19). 

Thus,  we  shall  approximate  B  in  (C19)  by  its  real  part  given  by 
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The  above  approximation  of  the  coefi icient  B  does  introduce  some  phase 
error  in  our  final  formula  (to  follow  shortly)  which  may  be  partially 
compensated  by  taking  the  absolute  value  of  A  in  (C18)  in  the  final 
formula.  This  compensation,  however,  is  not  seen  until  our  curve-fit 
formula  for  K+(a)  is  used  in  the  context  of  the  reflection  coefficient 
R(6.)  =  n/2n  K+(k)K+(-k  cose.)  . 

With  the  above  approximation  for  the  coefficients  A  and  B  duly 
noted,  we  may  write  from  (C17), 
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I A  |  [1  +  Br  (g)] 

K+°  (a) 

t  A 1  [  1  +  Br(g)+  C(±)] 


;  0  <  a  <  k 


;  -k  <  a  <  0 


(C22) 


which  leaves  only  the  coefficient  C  undetermined.  An  expression  for 
C  is  easily  obtained  by  matching  the  singularity  at  a  =  -k  in  (C22) 
with  the  actual  singularity  in  K(a)  at  a  =  -k  .  Using  the  relation¬ 
ship 

lim  K(a)  -  K+(k)  lim  K+(a)  (C23) 

a-*-  -k  a->  -k 


it  is  not  difficult  to  show  that 

1  -  | A | 2  (1  -  Bj) 

C  *  |a|4  (1  +  B  )  (C24) 

i  i  r ' 


Numerical  comparisons  between  our  approximate  expression  for  K+(a) 
in  (C22)  and  the  exact  formula  in  (Cll)  as  well  as  between  the  subsequent 
approximate  and  exact  reflection  coefficients  ( R ( 0 ^ )  =  n/2n  K+(k) 

K+(-k  cos6^))  are  found  in  Section  A  and  show  (C22)  to  be  a  very  good 
approximation  for  K+(a)  in  the  range  -k  <  a  <  k  for  ka  as  large 


as  1.0. 
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Appendix  D.  Limiting  forms  of  the  currents  on  cylindrical  antennas  near 
grazing  incidence  and  near  the  ends 


It  is  interesting  to  examine  the  limiting  forms  of  our  current 
expressions  for  semi-infinite  as  well  as  finite  length  cylindrical 
antennas  as  the  incident  angle  of  the  unifotm  plane  wave  approaches 
grazing  incidence  in  the  receiving  situation  or  as  an  end  is  approached 
in  either  the  receiving  or  transmitting  situations,  even  though  these 
limits  will  violate  the  basic  condition  placed  upon  our  theory  in  (8). 

It  is  noted  that  in  the  limits  to  follow,  the  value  of  the  electrical 
radius,  ka  ,  has  been  necessarily  taken  to  be  small.  This  condition  was 
imposed  to  avoid  complications  which  would  only  serve  to  obliviate  err 
original  purpose  of  seeking  a  qualitative  feeling  for  the  behavior  of  our 
formulas  beyond  their  expected  analytical  ranges. 


D.l  Semi- inf inite  receiving  antennc 

The  internal  current  on  a  semi-infinite  (0  <  z  <  -»)  cylindrical 
receiving  antenna  predicted  by  (43)  in  the  near-grazing  situation,  ■■  0 

may  be  expressed  as, 

{IS<»(®i;  z>  lint  ~  Eefi<h;  2>ei  1  as  0i  a-  0  ,  0  <  z  <  oo  (D.l) 

where. 
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The  internal  current  given  in  (D.l)  is  well  behaved  near  0^=0  and  goes 
to  the  physically  anticipated  limit  of  zero  at  0^  =  0  .  This  may  be 
attributed  to  the  fact  that  the  internal  current  expression  in  (43)  upon 


which  (D.l)  is  based  is  an  exact  expression  and  must  possess  the  correct 
physical  behavior. 


The  external  current  on  a  semi-infinite  (0  <  z  <  «•)  receiving 
antenna  is  given  by  the  combination  of  the  primary  receiving  current, 
I^(0^;z)  in  (15),  and  the  reflected  current,  Irefl^’2)  in  (25),  i.e., 

Igoo(6i;z)  =  Eq{V(0^;  z)  -  V(Oi;O)R(O.)U(0.;z) }  ;  0  <  6.  <  tt  , 

0  <  z  <  °°  (I).  3) 

Using  (9)  for  U(0^;z)  ,  (16)  for  V(0^;z)  and  (26)  with  (27)  for  R(0  )  , 
we  find  in  the  near-grazing  situation. 
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as  0^  -♦  0  , 

0  <  z  <  °°  (D.4) 


which  "blows  up"  at  0^  =  0.  Because  of  the  violation  of  the  condition 
in  (8),  this  behavior  should  not  be  unexpected.  However,  numerical  data 
in  this  instance  has  shown  that  the  misbehavior  is  limited  to  quite  small 
values  of  0^  ~  10°.  Similarly,  near  the  end  of  a  semi-infinite  receiving 
antenna  we  have  from  (D.3)  with  (9),  (16),  and  (26)  with  (27), 
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Br(l-cos0p  -  C(l-cos  0^) 
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8^(1  - cos0^)  -  C 
1  -  cos  0^ 


;  Tl/2  <  0  <  Tf 


and. 


2Cw  +  Itt  -  2  ?,n(sin  0  /2)  n 

Under  the  assumption  that  the  electrical  radius,  ka  ,  is  smai'  compared 
to  unity,  the  above  quantities,  F^e.)  and  F^O^  ,  will  also  be  quite 
small.  Hence,  the  current  in  (D.5)  will  be  at  a  relatively  small  level 
at  z  =  0  . 

Turning  to  the  other  form  of  U(e..;z)  in  (13),  we  have  for  the 
behavior  of  the  external  current  in  (D. 3)  using  (16)  and  (26), 
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;as  0i  •*  0  ,  0<z<°°  (D.8) 


which  goes  to  zero  at  0i  =  0  .  So  even  though  the  basic  condition  in  (8) 
is  violated  here,  the  correct  physically  expected  result  is  still  obtained. 
Near  the  end  of  the  semi-infinite  cylinder,  we  have  from  (13)  with  (16) 
and  (26)  in  (D.3), 
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;  as  z  -*■  0 


(D.9) 


where  F^(0^)  is  given  in  (D.6).  At  z  =  0  ,  the  external  current 
predicted  by  (D.9)  is  typically  a  very  small  value. 

D.2.  Semi-infinite  transmitting  antenna 

The  external  current  on  a  semi-infinite  (0  <  z  <  °°)  cylindrical 

transmitting  antenna  with  a  delta  function  voltage  source  at  z  =  zQ  is 

T 

given  by  the  combination  of  the  primary  current,  l^z^jz)  in  (19),  and 

T  ' 

the  reflected  current,  ^refl^O’2^  *n  *  **e*» 

I^z^z)  *  ll<*0**>  ~  I^(z0;0)R(tt)U(tt;z)  ;  0  <  z  <  °°  (D.10) 


The  behaviors  of  the  current  near  the  end  using  both  forms  of  U(tt;z)  in 
(9)  and  (13)  are  found  to  be, 
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10 

respectively.  At  the  end,  z  *  0  ,  the  currents  predicted  by  (D.ll)  and 
(D*12)  are  small,  the  latter  which  employed  (13)  for  U(tt;z)  being 
smaller  than  the  former  which  employed  (9)  for  U(ir;z)  . 

Although  small  values  of  current  were  obtained  at  the  end  in  both 
the  receiving  and  transmitting  cases  above  using  both  forms  of  IKO^z) 
in  (9)  and  (13),  none  of  these  results  agree  with  the  negative  of  the 
internal  penetrating  current  at  the  end,  discussed  in  Section  6,  which 
is  a  function  of  powers  of  ka  .  This  is  not  an  unexpected  result,  since 
our  theory  is  based  upon  the  condition  that, 

ft(z)  “  2  £n(2z/a)  >>  |&n[2kz  sin^(0^/2)|  ,  which  is  violated  m  these 
limits.  Nevertheless,  the  above  limiting  forms  are  for  the  most  part 
sufficiently  well-behaved  to  permit  useful  data  to  be  obtained  in  the 
finite  length  cylindrical  antenna  formulas  (discussed  in  Section  5)  over 
most  of  the  antenna  length  and  for  a  major  portion  of  the  possible  angles 
of  incidence. 


D. 3  Finite  length  receiving  antenna 

The  internal  current  near  the  z  -  -h  end  of  a  finite  length 
(-h  £  z  £  h)  cylindrical  receiving  antenna  predicted  by  (54)  with  (9) 
used  for  U(0i;z)  in  the  near-grazing,  Qi  ~  0°  ,  situation  behaves 
essentially  as. 
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Hence,  the  internal  current  near  the  z  3  -h  end  of  a  finite  length 
receiving  antenna  predicted  by  (54)  with  (9)  "blows  up"  at  grazing, 

8  »  0  ,  incidence  in  much  the  same  manner  as  the  external  current  given 
in  (D.4)  in  which  (9)  was  also  used  for  U(B^;z)  in  the  grazing  incidence 
limit.  A  similar  functional  behavior  to  that  in  (D.13)  is  obtained  in 
the  other  near-grazing  situation,  0^  ■+  it  ,  i.e.. 
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The  behaviors  of  the  internal  current  near  the  opposite  end  of  the  antenna, 
z  -h  ,  may  be  obtained  from  the  above  limits  with  appropriate  re¬ 
interpretation  of  the  incident  angle  of  the  uniform  plane  wave. 

The  limiting  behavior  of  the  external  current  distribution  on  a 
finite  length  receiving  antenna  given  in  (35)  with  (9)  used  for  U ( 0 ^ ; z) 
as  0.  ■*  0  is  given  by. 
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The  current  distribution  near-grazing  stated  above  shares  the  same 

functional  behavior  with  respect  to  0..  as  the  semi- inf inite  cylindrical 

antenna  limiting  formula  in  (D.4->,  i.e.,  it  "blows  up"  at  0i  =  0  . 

Although,  the  physically  expected  result  is  a  zero  current  at  grazing 

incidence,  0^  ■  0  ,  the  above  misbehavior  of  the  receiving  current 
formula  using  (9)  for  U(0^;z)  does  not  invalidate  this  formulation,  but 
can  be  shown  in  most  cases  to  limit  the  range  of  applicability  to  incident 
angles  in  the  range,  10°  <  0^  <  170°.  Because  of  the  symmetry  of  the 
receiving  formula  noted  in  (37),  a  similar  behavior  for  the  receiving 
current  is  obtained  in  the  limit  as  0^  *■  it  . 

Also  to  be  considered  is  the  limiting  form  for  the  finite  length 
cylindrical  antenna  current  at  the  ends.  Again  using  (9)  for 
in  (35)  we  obtain. 
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t 

where  F^(0^)  and  Fj^^)  are  given  in  (D.6)  and  (D.7)  and  are  small  for 
values  of  ka  that  are  small  compared  to  unity.  The  current  at  z  *  -h 
predicted  by  (D.19)  is  typically  a  small  value,  although  it  will  not 
correspond  to  the  negative  of  the  internal  penetrating  current  at  the  end 
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discussed  in  Section  6.  The  behavior  of  the  current  near  the  opposite 
end,  z  *  +h  ,  may  be  obtained  by  means  of  the  symmetry  relationship  in 
(37). 

Repeating  the  limiting  procedures  for  the  receiving  antenna  current 
in  (35)  using  (13)  for  U(0^;z)  ,  first  of  all  for  the  internal  current 
near  the  z  =  -h  end  in  the  near-grazing  situation,  yields. 
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where  f^(h;z)  is  given  in  (D.2).  Hence,  the  dominant  behavior  of  the 
internal  current  near  z  =  -h  on  the  finite  antenna  is  essentially 
given  by  the  behavior  of  the  current  near  the  end  of  a  semi-infinite 
antenna  near  grazing  as  given  in  (D.l).  A  zero  current  is  thus  obtained 
at  *  0  .  The  behavior  of  the  internal  current  near  z  *  -h  on  a 
finite  length  receiving  antenna  from  (35)  with  (13)  as  the  incident  angle, 
tends  to  the  other  grazing  angle,  0^  ■+  it  ,  is  given  by. 
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n  “  1 


The  behavior  of  the  external  current  on  a  finite  length  cylindrical 
receiving  antenna  predicted  by  (35)  using  (13)  for  U(6i;z)  is  given  by. 
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The  functional  behavior  of  the  current  in  (D.  23)  with  respect  to  the 
Incident  angle,  0^  ,  is  the  same  as  the  behavior  noted  in  (D.8)  which 
also  employed  (13)  for  U<©i; z)  .  Consequently,  the  physically  anticipated 
zero  current  is  obtained  in  (D.  23)  at  grazing  incidence,  8^  *  0  .  This 
must  be  judged  as  a  fortunate  happenstance,  since  our  theory  is  not 
expected  to  hold  for  very  small  incident  angles.  Also  by  means  of  the 
symmetry  relationship  in  (37),  a  similar  functional  behavior  can  be 
observed  as  the  opposite  grazing  incidence  angle,  0^  =  it  ,  is  approached. 

The  external  current  near  the  end  of  a  finite  length  cylindrical 
receiving  antenna  predicted  by  (35)  with  IKO^jz)  given  by  (13)  has  the 
limiting  form. 
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where  F^Q^)  is  given  in  (D.6).  (D.25)  typically  predicts  a  very  small 

value  at  z  *  -h  .  Although  it  will  not  correspond  to  the  negative  of 
the  internal  penetrating  current  at  the  end  discussed  in  Section  6.  The 
functional  behavior  of  the  current  near  the  opposite  end  at  z  =  +h  ,  can 
be  shown  to  be  similar  to  that  in  (D.25)  via  the  symmetry  relationship 
in  (40). 

Comparing  the  limiting  forms  of  the  currents  in  (D.  13) ,  (D.  15), 

(D.  17)  and  (D.  19)  with  those  in  (D.  20) ,  (D.  21) ,  (D.23)  and  (D.25)  corres¬ 
ponding  to  the  use  of  (9)  and  (13),  respectively,  for  U((L;z)  in  the 
finite  length  receiving  antenna  formula  in  (35),  the  latter  fOi.n  of 
U(0i;z)  must  be  judged  the  more  attractive  in  the  near  grazing  as  well 
as  the  end  vicinity  situations.  Not  only  is  the  correct  zero  current 
specified  at  grazing  incidence  when  (13)  is  used  for  ll(9^;z)  in  (35). 
but  the  current  at  the  ends  is  smaller  and  thus  more  manageable.  Both 
of  these  characteristics  are  particularly  important  in  the  case  of  the 
electrically  short  antenna  discussed  in  Section  8. 

D.4  Finite  length  transmitting  antenna 

Finally  we  consider  the  limiting  forms  of  the  external  current  at 
the  ends  of  a  finite  length  (-h  <_  z  ^  h)  cylindrical  transmitting 
antenna.  First  using  (9)  for  U(Tt;z)  in  (38)  we  find, 

IT(z0,z)  -  V0[U8(h+z0)  +  CT(h+z0)J 


— 

' 

Fx(tt)  -  [1  -F1(ir)J 

p  (tt)  .  12k<h 

2C  +in 

. 

L  w  J 

;  as  z  -*■  -h  (D.26) 


1L2 


where  F^(lt)  and  F2(10  are  given  in  (D.6)  and  (D.7)  and  are  normally 
small  for  values  at  ka  sufficiently  less  than  unity.  And  using  (13) 
for  U(tt;z)  in  (38)  yields  the  behavior, 

IT(z0;z)  ,  V0(Us(h  +  z0)  +  CT(h+z0)] 


[ 

r 

F1(ir)  -  [1-F10r)] 

- i2k(h  +  z) Hn [ 2k(h  +  z) ] 

2C  +i7i 
w 

1 

1 

;  as  z  -*■  -h  (D.27) 

Both  of  the  above  forms  of  the  transmitting  current  at  the  end,  z  *  -h  , 
are  typically  small,  the  second  one  using  (13)  for  U(tt;z)  being  the 
smaller.  Although  neither  limit  approaches  the  negative  of  the  internally 
penetrating  current  at  the  end  predicted  by  the  relationships 
discussed  in  Section  6.  A  functional  behavior  similar  to  the  above  for 
the  current  near  the  opposite  end  of  the  antenna,  z  =  +h  ,  may  be 
obtained  via  the  symmetry  relationship  in  (40) . 


-teL 


Ml 


L 


,l'x7V^ 
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Appendix  E.  Approximation  for  the  summation  in  the  expression  for  the 
current  at  the  cylinder  end 


Consider  the  summation  from  (43)  of  the  text  with  z  =  0, 
~  ik(Y0n  -  ik) 

S=  1  (Yn  +  i  k  cosTj  K+ ( iY0n  } 

n=l  On  On  i 


(El) 


Numerical  data  reveals  that  the  asymptotic  form  of  K+(a)  in  (46)  may  be 
used  for  K+(1Yq^)  with  little  error  provided  that  ka  <  1.  Then  (El) 


may  be  rewritten  as 

oo 

S  ~  ika  l 


[aY0n  "  ika] 


.  aYn  (aYn  +ika  cos  0.) 
n=l  On  On  i 


(E2) 


Since  we  are  considering  ka  <  1  and  since  p_  >  2.4048 

~  Un  — 

appropriate  Taylor  series  expansions  may  be  performed  on  the  terms  in  (E2) , 
leading  to 

S  ={(ika)  l  p  ^/2  -  (ika)2[i  +  cos  0  ]  £  p  5/2 

n-1  Un  1  n=l  un 


-(ika)3[|  -  ±  cos  6  -cos2  9  ]  l  p ~l'2 

n=l 


-(ika)^[-  ~  -  cos  et  +  |  cos‘  0i  +  cos3  0iJ  £  pQJ  2  +  •  •  •} 

n=l 


(E3) 


The  summations  in  (E3)  involving  the  zeroes  of  the  Bessel  function,  JQ, 


may  be  approximated  by 

oo  20  f 

V  -m/ 2  _  c  -m/2  .  .  .-m/2 

L  PQn  ~  4  pn«  +  J  I(x-1/4)tt]  dx  ; 


n-1 


n-1 


On 


20.5 


m  =  3,5,7, .. .  (E4) 
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The  integral  in  (E4)  approximates  the  original  summation  above  n  =  20  and 
utilizes  the  asymptotic  form  for  the  zeroes,  ~  (n-l/4)ir.  For  values 

of  ka  less  than  unity,  S  is  fairly  well  approximated  by  the  first  four 
summations  shown  explicitly  in  (E3).  Performing  these  summations  in  the 
approximate  manner  of  (E4)  allows  us  to  write  S  in  the  approximate  form, 

4  m 

S  =  y  A  (0 . )  (ika)  (E5) 

,  mi 
m=i 


where 

Al(9i)  =  °*5831  (E6) 

A2(0i)  =  -O.1364[|  +  cos(0.)]  (E7) 

AJ0J  =  -0.0498[-|-i  cos  0.  -  cos20  .  ]  (E8) 

j  i  o  z  i  i 

and 

A  (0J  *  -0.0198[~  cos  0.  +~  cos20.  +  cos30 .  ] 

4  i  lb  o  izi  i 

(E9) 


Appendix  F.  The  end  conductance  of  a  f i nite  length  cy  1  indr ica L  an tonna 


An  interesting  quantity  worth  mentioning  is  the  conductance  seen  by 

a  TMq^  mode  inside  and  near  one  of  the  ends  of  a  cylindrical  antenna. 

Employing  the  Wiener-Hopf  technique,  the  reflected  current  due  to  an 

1  ON  Z 

internal  current,  which  is  of  the  form,  e  ,  incident  upon  the  end 
of  a  semi-infinite  (0  <_  z  <_  £D)  cylinder  if;  approximately  given  by,  [31], 


ts°°  t  \  1 

Xrefl(z)  -  2 


k  +  iy 


ON 


ON 


K+“V  l 


K+(iW  z 


'On 


ON'  '0ny0n(Y0N+Y0n) 

n  =  1 


R(0ON)U(GON;z)  ;  0  <  z  <  *>  (FI) 


where  K  (a)  is  discussed  in  Appendix  C,  iY  is  the  propagation 
constant  of  the  TMQ^  circular  waveguide  mode  given  in  (45),  is  the 

ordered  zero  of  the  Bessel  function,  J  ,  and  finally  R(0,^.,)  and 

are  the  reflection  coefficient  and  canonical  integral  evaluated 
ON 

at  0ON  ”  cos  ^ (~iYojj/k)  given  in  (26)  and  (3),  respectively.  The 
first  term  in  (FI)  corresponds  to  internally  reflected  Tl^  circular 
waveguide  modes,  while  the  second  term  corresponds  to  the  externally 
transmitted  current  which  in  this  form  is  valid  under  the  condition  in  (8) 
From  (FI)  the  TM^  mode  current  reflection  coefficient  at  the  point, 
z  »  Si  ,  inside  the  cylinder  is  simply. 


ro> 


k  +  iy 


ON 


2y, 


ON 


V*W 


2  -2W 

e 


(F2) 


Through  a  transmission  line  analogy  [25],  we  may  write  the  end  admittance 
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where. 


YS>>  « 

e 


1 1 

1  '  r0N(£) 


ON 


(FI) 


'‘on 


4'rrk 

iYONn  : 


120  it 


(F4) 


is  the  analogous  transmission  line  characteristic  admittance  for  the 

TM_„  circular  waveguide  mode. 

UN 

By  instituting  the  multiple  reflection  concept  utilized  in  this 

report  (see  Section  5),  we  may  formulate  an  expression  for  the  external 

current  on  a  finite  length  (-h  _<  z  <_  h)  cylinder  from  the  initial 

externally  transmitted  current  in  (FI).  This  first  component  of  the 

external  current  as  well  as  all  of  the  higher  order  reflections  propa- 

i-kz 

gate  approximately  as  e  ,  so  that  the  corresponding  wave  incidence 
is  0^  ■  it.  The  constructed  external  current  distribution  on  a  finite 
length  cylinder  due  to  an  incident  internal  TM^  mode  current, 

Y0N(h+z) 

e  ,  is  given  by. 


Iext(z)  =  -R(0ON)U(0ON  ;h+z) 


R<eON)U(9ON;2h) 

1-  [ R( tt)  U(tt;  2h)  ]‘ 


R(tt)U(tt  ",  h  -  z) 


R(0qn)  U(0ON;  2h)  R(tt)  U(  tt;  2h) 
1  -  [R(TT)UU;2h)]2 


R(tt)  U(tt  ;  h  +  z)  ; 

-h  <  z  <  h 


(F5) 


With  the  knowledge  of  the  external  current  incident  upon  the  end, 
z  ■  -h  ,  from  (F5),  the  penetrating  current  near  this  end  may  be  formu¬ 
lated  from  the  transmission  characterization  in  Section  6  and  is 


given  by 
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■rext-KLnt.  x 
I  (z) 


R(eON)U(eON:2h)R(7,)U(,,;2h)  r 

—  I  T0n(,)e 


Y0n(h+z) 


1  -  [R(7r)U(n;2h)  ]' 


n  =  1 


z  >  -h 


(F6) 


where  and  iYQn  3r6  given  in  and  ^)  of  the  text.  The 

combination  of  the  penetrating  TMqn  mode  current  in  (F6)  and  the 
original  internally  reflected  TM^  mode  current  from  (FI)  with  respect 

Y0N^+Z^ 

to  the  incident  current,  e 
reflection  coefficient,  a  distance,  SL,  from  the  end,  z  =  -h,  given  by. 


yields  the  mode  cuirent 


r  rn  rs”m  k  R(0oN>u<0oN:2h)R(*>u^2h>  rfQ  '2W- 

r0Na)  r0Na)  +~n~  Y - 2 -  R(V’)e  (F7) 

N  n  Y0N  1  -  [R(TT)U(TT;2h)]^  ° 


The  equivalent  transmission  line  analogous  end  admittance  corresponding 
to  the  TM0^  mode  reflection  coefficient  in  (F7)  for  the  finite  length 
cylinder  is. 


YeU) 


^  +  1  ON^ 
1-r0N(X> 


ON 


(F8) 


where  YQN  is  given  in  (F4).  Figure  FI  shows  the  TM^  mode  end  conductance 
at  z  =  2a  -h,  i.e.,  £=2a,  for  a  finite-length  (-h  <  z  <  h)  cylindrical 
antenna  where  0(h)  =  2?_n(2h/a)  =10  as  a  function  of  the  electrical 
length  between  kh  =  0.05  and  kh  =  3.5  .  Also  shown  is  the  end 
conductance  of  a  semi- inf inite  cylinder  for  the  corresponding  electrical 
radii  between  ka  *  0.00067  and  ka  =  0.0472.  Note  that  the  end 
susceptance  at  this  depth  with  the  semi-infinite  as  well  as  the  finite 
length  cylinder  is  nominally  given  by  -lm{Y  f  from  (F4).  As  could  be 
expected  the  end  conductance  in  both  situations  i9  quite  small  and  the 
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Figure  FI.  End  conductance  for  the  TM  circular  waveguide  mode  at  a 
depth  equal  to  one  diameter,  i.e.,  2a,  inside  the  end  of 

a  tubular  cylinder. 

■  Effective  aperture  end  conductance 

from  eq.  (F9)  using  (54)  with  (13). 

-  —  Direct  (Wiener-Hopf )  end  conductance 

from  eq.  (F8)  with  (13)  used  for 

u(et;z). 

.  Semi- inf inite  cylinder  (electrical  radius. 

shown  in  top  scale,  corresponds  to  the 
electrical  radius  of  the  finite  length 
cylinder  for  £2(h)  =  10)  end  conductance 
from  the  effective  aperture  relation  in 
eq.  (F9)  using  (43)  or  the  direct  (Wiender- 
Hopf)  relation  in  eq.  (F3) ,  both  formulations 
yielding  similar  results. 


Note:  the  exact  value  of  the  term,  K+(iYQ1),  calculated  using  the 
formula  given  by  Mittra  and  Lee  [15,  Sec.  5-2, (3)]  has  been 
used  in  all  of  the  above  formulations.  Otherwise,  eq.  (27)  has 
been  used  to  calculate  K+(a)  for  -k  <_  a  k. 


END  CONDUCTANCE  (Mhos  xIO 
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finite  length  end  conductance  oscillates  above  and  below  the  semi-infinite 
end  conductance  beyond  the  first  resonant  length,  kh  -  ti/2  .  At  the 
smaller  values  of  kh  ,  the  finite  length  end  conductance  takes  an 
interesting  although  unlikely  diversion  from  the  expected  smoothly 
decreasing  behavior.  The  end  conductance  may  also  be  calculated  from  the 
formula  of  Chang  and  Rispin  [30]  which  is  based  upon  the  effective 
aperture  of  the  antenna  and  in  the  present  situation  may  be  written  as. 


G 

e 


sin  0 . 


d0. 

i 


(F9) 


sc  i 

| X0i (0i ; 2-)/eq  |  Is  the  magnitude  of  the  short  circuit  ™qj  mode  current 
at  l  and  is  equal  (if  we  neglect  the  multiple  reflections  in  the 
region,  -h  _<  z  <  -h+Jl)  to  twice  the  magnitude  of  the  internally 
penetrating  TM^  mode  current  given  in  (43)  for  the  semi-infinite 
cylinder  and  in  (54)  for  the  finite  length  cylinder.  In  order  to  neglect 
the  perturbation  of  the  external  current  and  the  multiple  reflections 
of  the  TMqj  mode  within  the  end  of  the  cylinder  due  to  the  short  circuit 
at  z  *  £-h  it  is  necessary  to  require  that  £  >  2a  in  this  conductance 
formulation.  The  semi-infinite  and  finite  length  end  conductances  obtained 
from  (F9)  for  an  antenna  where  0(h)  =  2Jtn(2h/a)  =  10  is  also  included 
in  Figure  FI.  The  semi-infinite  cylinder  end  conductances  from  both  the 
direct  and  effective  aperture  methods  in  (F3)  and  (F9)  with  (54),  respec¬ 
tively,  are  for  all  practical  purposes,  the  same.  And  the  agreement 
between  the  two  finite  length  end  conductances  in  Figure  FI  at  the  larger 
values  of  kh  is  very  good.  However,  at  the  smaller  values  of  kh  ,  the 
effective  aperture  end  conductance  from  (F9)  with  (54)  continues  to 
decrease  smoothly  for  shorter  electrical  lengths  as  opposed  to  the  odd 
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behavior  of  the  directly  formulated  end  conductance  in  (F8) .  We  accept 
the  validity  of  the  effective  aperture  result  on  the  grounds  that  it 
depends  only  upon  the  magnitude  of  the  associated  current  which  is 
assumed  to  be  fairly  accurate,  since  the  magnitude  of  the  receiving 
current  on  the  electrically  short  cylindrical  antenn-'  pr-di  te bv  our 
formulas  compared  very  well  with  an  accepted  existing  theory  in 
Section  8.  On  the  other  hand,  the  direct  formulation  in  (F8)  depends  on 
the  phase  of  the  reflected  current,  which  we  know  is  not  accurately 
determined  by  our  formulas  applied  to  the  electrically  short  cylindrical 


antenna . 
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Appendix  G.  Electrically  thick  antennas  on  the  order  of  a  few  wavelengths 
in  length 

The  application  of  our  theory  to  the  electrically  thick  (up  to  ka  =  1) 
cylindrical  antenna  whose  length  is  on  the  order  of  a  few  free  space 
wavelengths  is  based  upon  observations  in  Section  4  which  imply  that  the 
basic  condition  in  (8)  may  be  relaxed  for  the  larger  values  of  ka  to 
the  point  where  (l(z)  =  2  £n(2z/a)  need  only  be  somewhat  larger  than 
|lln[2kz  sin^(0^/2)]|  to  achieve  reasonably  accurate  results  for  the 
primary  and  secondary  currents  on  the  antenna.  Although  the  application 
of  our  theory  (using  (13)  for  U(9^;z)  )  has  already  been  demonstrated 
for  the  electrically  thick  antenna  in  Section  7,  we  shall  present  results 
in  this  appendix  for  similar  and  even  shorter  antennas,  based  upon  our 
formulas  using  (9)  for  U(0i;z)  .  And  again  we  note,  that  for  all  of  the 
antennas  to  which  our  theory  is  applicable,  the  total  (external  and 
Internal)  current  distribution  is,  for  all  practical  purposes,  given  by 
the  external  current  distribution  alone,  except  in  the  near  vicinity  of 
the  ends, 

G.l  Electrically  thick  receiving  antenna 

The  current  distribution  on  an  electrically  thick,  ka  *  1  ,  finite 
length,  kh  *  3ir  ,  cylindrical  antenna  as  calculated  from  (35)  using  (9) 
for  UCS^jz)  is  shown  in  Figure  G1  for  the  incident  angles  (of  the 
uniform  plane  wave)  of  0i  -  tt/6  ,  it/ 3  ,  and  ti/2  .  Normal  incidence 
data  given  by  Wu,  et  al.  [35]  based  upon  the  theory  of  Kao  [18]  for  the 
azimuthally  uniform  z-directed  current  is  also  shown  in  this  figure 
for  comparison.  And  the  resulting  agreement  between  the  two  sets  of 
normal  incidence  (0^  ■  ir/2)  data  is  excellent  (except,  of  course,  near 
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Che  ends)  and  better  than  the  agreement  which  was  previously  obtained  and 
showry  in  Figure  21  in  which  the  other  approximate  form  of  ll(8.;z)  in  (13) 
was  used  in  the  receiving  current  expression  given  in  (35).  And  even 
further.  Figure  G2  shows  the  receiving  current  distribution,  as  calculated 
from  (35)  using  (9)  for  U(0^;z)  ,  on  an  electrically  thick,  ka  =  1  , 
antenna  whose  electrical  length  is  one  half  that  of  the  antenna  in  the 
previous  figure,  i.e.,  kh  =  3rr/2  .  The  agreement  between  our  normal 
incidence  (0^  “  tt/2)  result  and  the  numerically  obtained  normal  incidence 
data  of  Kao  [18]  is  again  excellent.  An  attempt  to  apply  our  theory  to 
this  antenna,  (ka  ■  1  and  kh  =  3tt/2  )  ,  using  the  other  approximate 
formula  for  U(0^;z)  ,  which  is  given  by  (13),  met  with  much  less  success. 
However,  of  all  the  many  examples  studied  in  the  preparation  of  this  report, 
it  was  only  for  these  electrically  thick  antennas  did  the  use  of  (9)  for 
U(0^;z)  show  a  decided  advantage  over  the  use  of  (13)  for  U(0i;z)  in 
our  finite  length  antenna  formulas.  Hence,  our  theory  involving  the  use 
of  (13)  for  U(0^;z)  was  featured  in  the  text  of  this  report.  And 
although  a  near-grazing  incidence,  0^  =  tt/36  ,  distribution  was  included 
in  the  receiving  antenna  examples  in  the  text,  no  near  grazing  results 
have  been  shown  in  Figures  Gl  and  G2 ,  due  to  the  misbehavior  of  our 
receiving  current  formula  for  small  incident  angles  (see  Appendix  D)  when 
(9)  is  used  for  U(6^;z)  . 

G. 2  Electrically  thick  transmitting  antenna 

A  recalculation  of  the  input  conductance  for  the  tubular  antennas 
considered  in  Figure  29  (which  used  (13)  for  U(ir;z)  in  (41)  for  the 
calculation),  this  time  using  (9)  for  U(7r;z)  in  (41)  results  in  the 
curves  shown  in  Figure  G3.  Data  from  the  numerical  approach  of  Chang  [27] 
and  [28]  and  the  experimental  method  of  Hartig  [29]  which  were  used  for 
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Figure  G2.  Current  distributions  on  a  three-halves  wavelength  long, 

electrically  thick,  cylindrical  receiving  antenna  normalized 
to  the  incident  electric  field  and  the  wavelength,  i.e., 

IR(0i;2)/XeJ  . 


from  eq.  (35)  with  (9) 
used  for  U(6  ;z) 

Numerically  evaluated  data 
from  Kao  [ 18] . 
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comparison  in  Figure  29  is  also  included  in  Figure  G3.  A  close  inspection 
and  comparison  of  Figures  29  and  G3  reveals  about  the  same  leveL  oi  agreement 
between  our  theory  and  the  others  in  the  latter  figure  where  (9)  was  used 
for  U(tc;z)  .  However,  as  explained  in  the  previous  sub-section,  (13) 
was  used  for  U(tt;z)  in  the  text  examples,  since  it  had  shown  a  wider 
range  of  applicability. 

Further,  in  Figure  G4  we  have  repeated  the  conductance  data  from 
Figure  30  for  a  center-fed  tubular  antenna  where  h/a  =  10  and  included 
results  obtained  from  (41)  in  which  (9)  has  been  used  for  U(ir;z)  .  We 
note  that  in  this  figure  the  two  forms  of  U(7T;z)  (given  by  (9)  and  (13)) 
when  used  in  the  admittance  formula  given  in  (41)  yield  slightly  different 
results,  although  both  conductance  curves  are  consistent  with  the  numerical 
data  of  Harrington  and  Mautz  [26] .  Figure  G5  is  identical  with  Figure  31 
which  showed  the  input  conductance  for  an  off-center  fed,  z„  =±h/2  , 

ill  o 

tubular  antenna  where  h/a  *  10  ,  except  for  the  additional  conductance 
data  obtained  from  (41)  this  time  using  (9)  for  U(tt;z)  .  And  again  all 
three  sets  of  data  are  consistent. 
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Figure  G4.  Input  conductance  of  a  center  fed  cylindrical  antenna  as  a 
function  of  the  normalized  length,  2h/A  . 
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Figure  G5,  Input  conductance  of  an  off-center  fed  cylindrical  antenna 
as  a  function  of  the  normalized  length,  2h/A  . 
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Appendix  H .  Far  field  radiation  from  a  cy lindrical  _tr ansmitting  ant enna 
The  far-field  radiation  pattern  of  a  finite  length,  -h  <'  z  h  , 
cylindrical  antenna  having  a  delta  function  source  of  strength,  volts, 

at  z  =  Zq  (similar  to  the  antenna  shown  in  Figure  16c)  can  be  written  in 
terms  of  [33,  Sec.  2-10], 


Eg(r,0,4>)  =  — |^sin0  — 


ikr 


-h 


—  7T 


,1.  ,  ,,  -ikz’cosO 

,z  )e 


-ikasinOcos  <$>' 
e 


ad <p  '  dz  '  ; 


kr  ■->  1 


(H.  1) 


T 

where  Jz(a,<f>,z)  is  the  z-directed  component  of  current  density  on  the 

cylinder.  E.(r,0,rf))  is  the  8-directed  (  0  is  measured  from  the  positive 
0 

z  axis  in  Figure  16c)  electric  field  at  the  far  field  point  (r,u,<p)  i;.  a 

spherical  coordinate  system  coincident  with  the  implied  cylindrical  coor- 

T 

dinate  system  of  Figure  16c.  Since  Jz(a.<p,z)  is  uniform  about  the  cylin- 

T  T 

der,  we  take  the  total  current  as  I  (z^;z)  =  27TaJz(a,<j>,z)  and  performing 
the  <J)'  integration,  we  get. 


Eg(r,0,(}))  *  — ^  sin  0  Jg(ka  sin  0) 


ikr 


tT,  , ,  -ikz  'cost1  ,  , 
I  (zQ;z')e  dz 


kr  >>  1 


(H.2) 


where  Jq  is  the  Bessel  function  of  the  first  kind.  Although  our  approxi- 

T 

mate  expression  for  I  (Zgjz)  in  (38)  may  be  used  in  (H.2)  to  determine 
the  far  field,  the  integration  of  some  of  the  terms  requires  further 
approximations  and  leads  to  a  more  complicated  result  than  is  actually 
necessary. 

The  far  field  expression  in  (H.2)  may  be  stated  in  a  more  convenient 
manner  if  we  first  consider  a  similar  receiving  antenna.  The  current  dis¬ 
tribution  on  an  antenna  used  as  a  receiving  element  can  be  found  by 
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integrating  the  response  of  the  antenna  to  a  unit  voltage  impulse,  i.e., 

I T(z';z)/Vq,  with  the  differential  voltage. 

pi/  o  i\j  •  ■  a  i  pi  ikz 'cos(ti-0)  , 

Eo(a,0,z  )dz  =  sm  8  J.(ka  sin0)Eoe  dz  (H.3) 

o  0  o 


over  the  range,  -h  _<  z  <  h  .  Here  Eg(a,0,z’)dz'  corresponds  to  the 
azimuthally  uniform  differential  component  of  voltage  applied  to  the 
antenna  when  illuminated  by  a  uniform  plane  wave  having  a  O-directed  com¬ 
ponent  of  electric  field.  Eg  ,  and  incident  at  an  angle,  n  -  0  ,  with 
respect  to  the  positive  z-axis.  Noting  that  the  exact  form  of  1  (z';z) 

.  *p  *p 

would  satisfy  the  reciprocal  relationship.  I  (z';z)  =  I  (z;z')  ,  [34. 

Sec.  9-10],  the  receiving  current  distribution  can  be.  written  as. 


E1  r  h  , 

IR(n-0  ;  z)  =-sin  0  Jn(ka  sinQ)  ~  lT(z;  z '  )e“lkz' cos0  dz’  (H.4) 

0  J-h 


Using  (H.4)  with  z  =  z^  in  (H.2)  allows  the  far  field  to  be  written  as. 


. ,  ikr  I  (tt  —  0  ,  z  ) 

E0(r.0,^irVVO - “1“^  kr>>1 

Ee 


(H.5) 
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